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ABSTRACT 

We use interferometric ^^CO(l-O) observations to compare and contrast the extent, surface 
brightness profiles and kinematics of the molecular gas in CO-rich ATLAS'^'^ early-type galax- 
ies (ETGs) and spiral galaxies. We find that the molecular gas extent is smaller in absolute 
terms in ETGs than in late-type galaxies, but that the size distributions are similar once scaled 
by the galaxies optical/stellar characteristic scalelengths. Amongst ETGs, we find that the ex- 
tent of the gas is independent of its kinematic misalignment (with respect to the stars), but 
does depend on environment, with Virgo cluster ETGs having less extended molecular gas 
reservoirs, further emphasizing that cluster ETGs follow different evolutionary pathways from 
those in the field. Approximately half of ETGs have molecular gas surface brightness profiles 
that follow the stellar light profile. These systems often have relaxed gas out to large radii, 
suggesting they are unlikely to have had recent merger/accretion events. A third of the sam- 
ple galaxies show molecular gas surface brightness profiles that fall off slower than the light, 
and sometimes show a truncation. These galaxies often have a low mass, and either have dis- 
turbed molecular gas, or are in the Virgo cluster, suggesting that recent mergers, ram pressure 
stripping and/or the presence of hot gas can compress/truncate the gas. The remaining galax- 
ies have rings, or composite profiles, that we argue can be caused by the effects of bars. We 
investigated the kinematics of the molecular gas using position-velocity diagrams, and com- 
pared the observed kinematics with dynamical model predictions, and the observed stellar and 
ionised gas velocities. We confirm that the molecular gas reaches beyond the turnover of the 
circular velocity curve in ?s70% of our CO-rich ATLAS'^'^ ETGs, validating previous work 
on the CO TuUy-Fisher relation. In general we find that in most galaxies the molecular gas is 
dynamically cold, and the observed CO rotation matches well model predictions of the circular 
velocity. In the galaxies with the largest molecular masses, dust obscuration and/or population 
gradients can cause model predictions of the circular velocity to disagree with observations of 
the molecular gas rotation, however these effects are confined to the most star-forming sys- 
tems. Bars and non-equilibrium conditions can also make the gas deviate from circular orbits. 
In both these cases one expects the models circular velocity to be higher than the observed 
CO velocity, in agreement with our observations. Molecular gas is a better direct tracer of the 
circular velocity than the ionised gas, justifying its use as a kinematic tracer for Tully-Fisher 
and similar analyses. 

Key words: galaxies: elliptical and lenticular, cD - galaxies: evolution - galaxies: ISM - 
ISM: kinematics and dynamics - ISM: molecules - ISM: evolution 



1 Introduction 

The buildup of galaxies through hierarchical merging has 
emerged in recent years as the predominant theory of galaxy evolu- 
tion. In this paradigm, the collapse of dark matter haloes proceeds 
in a "bottom-up" fashion, smaller structures forming first and merg- 
ing to create larger systems (e.g. [Springel et al.|[2005^ . It seems 
however that galaxy star-formation histories do not share the same 
"bottom-up" evolution. Massive early-type galaxies (ETGs) are the 
end-point of hierarchical galaxy evolution; they are normally con- 
sidered to be 'red and dead' , without significant star-formation, and 
containing mainly old stellar populations (e.g. [Bower, Lucey & El-| 
|lis|1992[ l. Lower mass spiral and irregular galaxies, which generally 
exist in less dense environments and have undergone fewer merg- 
ers, are usually gas rich and actively star-forming (e.g. [Roberts &| 
[Haynes[1994| . 

Some ETGs do, however, have molecular (and atomic; see 
[Morganti et al.[2006[|Oo"sterloo et al.[2"oTo|[Serra et al.[2011[ Paper 
XIII) gas reservoirs. Molecular gas was first detected in ETGs by 
[Wiklind & Rydbeck[ p986) and [Phillips eTaT] (1987) 1. Since then 
various surveys have detected molecular gas in ETGs that are in- 
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frared bright (e.g. [Wiklind & Henkel|1989i[Sage & Wrobel|1989[ 



[Knapp & R upen'1996 1, or that have dust obscuration in optical im- 
ages (e.g. Wang, Kenney & Ishizuki[1992^ . It was not until recently 
however that unbiased samples of ETGs were searched for molec- 
ular gas (e.g. [Welch & Sage"2003', 'Sage, Welch & Young[[2007l 
[Combes, Young & Bu reau 2007 1. Young et al. ( 201 1] hearafter Pa- 
per IV) conducted the largest of these studies to da te, on the com- 
plete, volume-limited ATLAS'^° sample of ETGs (Cappellari et al. 
[20TT] hereafter Paper I). This work found that 22% of ETGs in 
the local volume (out to 42 Mpc) have a detectable molecular gas 
reservoir (Mh2 ^ 1 x lO^^M©). It is thus important to understand the 
origin, properties and fate of this gas, to understand downsizing and 
the evolution of massive galaxies. 

To explain the tight correlations between the properties of 
ETGs (e.g the fundamental plane, [Dressier et al. 1987; Djorgovski] 
[& Davis[1987[ the colour-magnitude relation, Baldry et al. 2004 1, it 
is thought that the cold inter-stellar medium (ISM) must be removed 
or destroyed (or potentially made stable; e.g. [Martig et al.[[2009t , 
such that star-formation is abruptly quenched. How one reconciles 
the requirement for quenching with the observation of molecular 
and atomic gas reservoirs in some ETGs is an open question. In 
Paper X of this series ( [Davis et al.][20TT) , we showed that much 
of the molecular gas in field ETGs is kinematically misaligned and 
likely has an external origin. Thus, consistent with the quenching 
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scenario, these galaxies could have been red and dead in the past, 
and have recently acquired more gas. Virgo cluster ETGs also have 
molecular gas however, without replenishing this externally. 

[Khochfar et al.| \2Ql \\ (hereafter Paper VIII) suggest, from a 
modelling perspective, that star formation and hence molecular gas 
are important in the transformation of galaxies already on the red- 
sequence (e.g. in turning slow rotators into fast rotators). In order 
to test such a hypothesis one needs to understand the state of the 
molecular gas, and its sources. 

The molecular gas in ETGs exists in different conditions than 
those within the discs of spiral galaxies. For example, ETGs have 
different gravitational potentials (affecting e.g. the |Toomre||l964| 
stability parameter Q), different elemental abundances (high metal- 
licity and a-element over-abundance), higher prevalence of hot gas 
(10® K) and active galactic nuclei (AGN). Old stellar populations 
also create a harder inter-stellar radiation field (e.g. the UV-uptum 
phenomenon; see Yi 2008: B ureau et al.||20l7) . Furthermore, as 
mentioned above, many galaxies have kinematically misaligned gas 
with an external origin, rotating along a different axis than the ma- 
jority of the stars. Spiral galaxies have had their star-forming gas 
in situ for many gigayears (even if it is also being topped up ex- 
ternally; e.g. |Fratemali & Binney|2008^ . All these differences may 
affect the properties of the gas, including its extent, ability to form 
stars, and kinematics. 



To investigate these issues, it is vital to have spatially-resolved 
data on the distribution and kinematics of the molecular gas. Until 
recently, only a small number of ETGs had been mapped with mm- 



interferometers (e.g.|Wrobel & Kenney 1992; Young 2002 


|Young, 


Bureau & Cappellari||2008i ,Schinnerer & Scoville||2002l 


Crocker 


et aI.|2008||2009||Wei et al.|20I0||Crocker et al.|201 IJ, hence draw- 



ing conclusions based on a few objects has been difficult. In this 
paper, we attempt to understand the similarities and differences be- 
tween the extent and distribution of the molecular gas in spirals 
and ETGs. For this we will use a sample of 41 ETGs with molec- 
ular gas mapped as part of the ATLAS'^° project (Alatalo et al., in 
prep; hereafter A 12). Furthermore, we continue a theme presented 
in Paper V of this series ( [Davis et al.|2011^ , evaluating the use of 
molecular gas as a circular velocity tracer in ETGs, this time in a 
spatially-resolved manner. 

In Section 2, we present our galaxy sample and outline the ob- 
servations and models our work is based on. In Section 3 we present 
a study of the extent of the molecular gas in our ETGs, and compare 
this to the extent of the gas in spirals, and this extent varies with 
gas origin and environment. In Section 4 we examine the surface 
brightness profiles of the molecular gas in ETGs. In Section 5 we 
evaluate the use of CO as a tracer of the circular velocity in ETGs, 
and how the observed CO rotation compares with circular veloc- 
ity predictions from mass models and observed stellar and ionised 
gas rotation velocities. We discuss these results in Section 6, before 
concluding in Section 7. 

2 Data and Methods 
2.1 Sample 

We present results based on a sample of 41 ETGs that have ^^CO 
aperture synthesis imaging available, either from the ATLAS^° 
project (A12) or the literature. As part of the ATLAS'^° survey, all 
CO(l-O) detections from Paper IV with an integrated flux greater 
than 19 Jy km s~^ that do not have interferometric data available in 
the literature have been observed with the Combined Array for Re- 
search in mm-wave Astronomy (CARMA; |Bock et al.|2006| l. Full 



details of this interferometric survey can be found in AI2, but we 
summarise the observations in Section|2j2] 

A total of 28 galaxies included in this work have been 
mapped with CARMA (see Table We also include galaxies for 
which data are already available from the literature, mostly from 
the SAURON survey l |de Zeeuw erS^[2002) follow-ups. These 
are NGC0524, NGC2685, NGC2768, NGC3032, NGC3489, 
NGC4150, NGC4459, NGC4476, NGC4477 NGC4526, and 
NGC4550. Full references for the observations used are listed in 
Table [T] Furthermore we include NGC2697 and NGC4292, ob- 
served with CARMA as part of this project but later removed from 
the final ATLAS^'^ sample (due to not making the distance and/or 
magnitude cut). These two objects are nevertheless classified as 
ETGs, and by including them in the sample we are unlikely to in- 
troduce any significant bias. In total, this results in a sample of 41 
galaxies that have CO(l-O) interferometric data. 

This sample of mapped galaxies includes the brightest ~2/3 
of the CO detections from Paper IV. The mapped galaxies are thus 
biased towards higher molecular gas masses, but are statistically 
indistinguishable (using a Kolmogorov-Smirnov test) from the full 
sample of CO-detected systems in terms of their gas fractions (the 
ratio of molecular mass to stellar mass), group scale environment, 
and host galaxy properties e.g. velocity dispersion, luminosity, ef- 
fective radius (Paper I) and specific angular momentum parameter 
(Ai;; |Emsellem et al.|201 1[ hearafter Paper III). The mapped galax- 
ies also follow the same CO TuUy-Fisher relation (Paper V). Other 
than molecular gas mass, the only parameter in which our sample 
may have some bias is cluster membership. Because the Virgo clus- 
ter is nearby, its members are relatively bright in CO(I-O), and we 
have mapped a greater percentage of cluster members than field 
galaxies. However, we do not expect this to affect our conclusions. 
Overall, we thus consider that our mapped CO-rich ETG sample 
is reasonably free of biases, although it is worth remembering that 
we do not fully sample the parameter space at small molecular gas 
masses. 

2.2 CARMA data 

Observations of the 28 sample galaxies were taken in the D-array 
configuration at CARMA between 2008 and 2010, providing a spa- 
tial resolution of 4-5". Full details of this interferometric survey can 
be found in AI2 but we summarize the observations here. CO(l-O) 
was observed using a narrow-band correlator configuration, provid- 
ing at least 3 raw channels per 10 km s^^ binned channel whilst 
ensuring adequate velocity coverage for all galaxies. Bright quasars 
were used to calibrate the antenna-based gains and for passband cal- 
ibration. The data were calibrated and imaged using the 'Multichan- 
nel Image Reconstruction, Image Analysis and Display' (MIRIAD) 
software package jSault, Teuben & Wright 1995 1. Comparison of 
the CARMA total fluxes with the IRAM 30m single-dish obser- 
vations from Paper IV shows that within the measurement errors 
we always recover all the flux in our CARMA maps, and often 
have higher fluxes due to detecting gas outside the 30m telescope 
beam (A12). The sensitivity of our observations vary from 3-30 
mJy/beam (dependant on the channel size used in the final cube; 
average lOkms"^). This gives us a 3(t surface density limit of be- 
tween 3 and 30 Mq pc~^ with a median of 14 Mq pc^^ (if one 
uses a standard galactic Xco factor as in Paper IV). 

The CO integrated intensity maps we use in this paper are pre- 
sented in AI2. They were produced by the masking method. The 
cleaned image cube was smoothed spatially with a gaussian 1.4 
times larger in diameter than the beam, and in velocity with a ban- 
ning width of 3 channels, and then this smoothed cube was clipped 
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at ~2.5 times the rms noise in a channel. The chpped cube was then 
used as a mask to define a three-dimensional volume in the original, 
unsmoothed cleaned cube, over which we took the zeroth moment 
(see |Regan et al.|2001^ . For some of the analysis in this paper, un- 
dipped moment maps were created from the reduced data cubes of 
A 12, to avoid biasing against faint emission. 

2.3 SAURON IFS data and Jeans Models 

The SAURON observations of the ATLAS^° galaxies and the ex- 
traction of the stellar kinematics are described in detail in Paper 
I. In brief, for each target, individual datacubes were merged and 
analysed as described in |Emsellem et al., ( 2004 l, ensuring a mini- 
mum signal-to-noise ratio of 40 per spatial and spectral pixel using 
the binning scheme developed by |Cappellari & Copin|p003) . The 
SAURON stellar kinematics were derived using a penalized pixel 
fitting routine ( [Cappellari & Emsellem|2004^ , providing paramet- 
ric estimates of the line-of-sight velocity distribution for each bin. 
During the extraction of the stellar kinematics, the GANDALF code 
( |Sarzi et al.|2006l l was used to simultaneously extract the ionised 
gas line fluxes and kinematics. 

In this paper we utilize axisymmetric leans anisotropic dynam- 
ical modeling (JAM; .Cappellari. 2008 j of the ATLAS^°galaxies. 
Some examples of using this modeling approach with SAURON 
integral-field kinematics are presented in jScott et al.| ( |2009l). For the 



ATLAS''^ survey a multi-Gaussian expansion (MGE; 



Emsellem. 



Monnet & Bacon 19941 was fitted to the SDSS (Abazajian et al. 



,2009 ) or Isaac Newton telescope photometry (Scott et al., in prepa- 
ration). The MGEs were then used to construct JAM models for all 
260 ATLAS^° galaxies (see [Cappellari et aL|[2oTo |, fiting to the 
SAURON stellar kinematics (Paper I). The models have three free 
parameters, the inclination (ijam), the mass-to-light ratio (M/L) 
assumed to be spatially constant, and the anisotropy also assumed 
to be spatially constant. 

From each mass model, we have calculated the predicted cir- 
cular velocity curve in the equatorial plane of the galaxy. It is worth 
bearing in mind that where the CO is misaligned from the plane 
of the galaxy, using the circular velocity curve from the equatorial 
plane would only be valid if the potential was spherically symmet- 
ric. As these models do not include dark matter, the circular velocity 
often declines at large radii, and one must be careful to compare cir- 
cular velocities only where the models are constrained by the IFS 
data. Even if substantial amounts of dark matter are present in the 
central parts of these galaxies, where the models are constrained by 
the IFU data the circular velocity predicted will be correct ( jCappel-| 
|lari et al.,2006 ). As in this paper we only consider the inner-parts of 
the galaxies, where the models are fully constrained, these consid- 
erations are unlikely to affect our results. 

3 Molecular Gas Extent 

In this section we attempt to understand if the radial extent 
of the molecular gas in ETGs is different from that in spirals. One 
can imagine mechanisms that would make the gas in ETGs less 
extended than the gas in spirals. External gas, for instance, could 
fall into the galaxy centre, and create inner discs and rings ( [Eliche-| 
[Moral et al.||201l) . Furthermore, the interaction of kinematically- 
misaligned gas with stellar mass loss from stars could cause the gas 
to loose angular momentum and become less extended over time. To 
look for signatures of these processes, we measured the extent of the 
gas in 41 CO-mapped ATLAS^"^ ETGs, and a comparison sample 
of 44 spirals from the Berkeley Illinois Maryland Array Survey of 
Nearby Galaxies (BIMA-SONG; |Helfer et al.|2003^ . 



3.1 Redshifting BIMA-SONG 

BIMA-SONG mapped CO(l-O) molecular emission in a volume 
(and magnitude) limited sample of 44 nearby spiral galaxies. The 
survey included every galaxy of Hubble type Sa - Sd (except M33 
and M31), with declinations 5 > -20°, visual magnitudes B < 1 1.0, 
velocities Vhd < 2000 km s^^, and inclinations i < 70°. The sur- 
face brightnesses sensitivity of the BIMA-SONG survey is very 
similar to our ATLAS^° survey (with a 3a limit of «20 Mq pc"^ 
when adjusted to our preferred Xco; [Heifer et al.|2003^ . As such 
it is an ideal survey to compare with our ATLAS'^^galaxies, which 
are also drawn from a complete volume limited survey. We chose 
not to use the HERACLES survey ( |Leroy et al.|2009| for our com- 
parison, because the sample was not selected in a volume limited 
manner, the angular resolution achieved is 2-3 times worse than our 
observations (13") and the galaxy sample is small (18 objects). We 
do however compare qualitatively with the HERACLES results in 
Section leXn 

The molecular emission in these objects was mapped with a 
typical spatial resolution of 6" (360 pc at 12 Mpc, the average dis- 
tance of the sample galaxies; |Helfer et al.|2003[ l. The array configu- 
rations used had a minimum baseline of «8m, meaning the largest 
detectable scales are similar to our CARMA data (which was ob- 
served with a minimum baseline of 11m). We caution, however, 
that |HeIfer et al.|f2003^ showed that they often resolved out signif- 
icant fractions of the flux in these spiral galaxies. In order to min- 
imise the effect of this missing flux we used the data cubes which 
included zero-spacing information (from on-the-fly mapping with 
the NRAO 12m telescope), where available (24/44 objects). In the 
other galaxies we were unable to correct for the distribution of this 
missing flux, and thus caution that the sizes derived here should 
be strictly considered as lower limits. If we only consider galaxies 
with zero-spacing information available, the trends reported below 
do not change however. 

Another difficulty we faced when comparing the samples was 
that the ATLAS^'^ sample galaxies are on average twice as dis- 
tant (24.5 Mpc) as the BIMA-SONG galaxies. To overcome this 
problem we artificially 'redshifted' the BIMA-SONG galaxies to 
the average distance of the ATLAS''^ galaxies, using the procedure 
outlined below. 

Firstly, we obtained the clean cubes from the BIMA-SONG 
archive. The angular size of the beam when observing a distant 
galaxy appears larger relative to the galaxy size, so each BIMA- 
SONG cube was convolved with a Gaussian to create the correct 
relative beam size. The size of the major and minor axes of the con- 
volving Gaussian required in order to reach the correct final beam 
were calculated by assuming that the beams add in quadrature: 



B° 



d 



24.5 Mpc 



(1) 



where -Bmfn.maj is the original beam size, -Bmin.maj the convolving 
beam size, d is the distance to the BIMA-SONG galaxy and 24.5 
Mpc is the average distance to the ATLAS-3D galaxies. 

Secondly the clean BIMA-SONG cube was scaled to the flux 
expected at the larger distance. This part of the procedure is unlikely 
to effect our ability to detect these objects, as surface brightness 
sensitivity is independent of distance. 

Thirdly, a noise cube was created with the Miriad task IM- 
GEN. It was then scaled to the correct level, so that when added to 
the BIMA-SONG clean cube the resultant cube had the same RMS 
surface brightness threshold as a typical CARMA observation of 
the ATLAS^° galaxies (15 Mq pc~^ in a 10 km s~^ chan). If the 
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BIMA-SONG cube already has an RMS greater than this level, no 
additional noise was added. This only applies in some cases, how- 
ever, and we do not expect this to bias our results. 

Finally, a moment zero map of the "redshifted" cube was cre- 
ated, using the same procedure outlined in Section [T2| The pixel 
size was set so the map is Nyquist sampled (3 pixels across the ge- 
ometric mean of the beam size). 

3.2 Measuring gas extents 

The molecular gas extent of each galaxy was estimated from its 
ATLAS^° or redshifted BIMA-SONG moment zero map. The di- 
ameter was measured in the direction of the largest molecular gas 
extent, a level of approximately three times the RMS noise in the 
moment zero maps. As the 3cr surface brightness limits are approx- 
imately the same in the two surveys we should be probing simi- 
lar fractions of the gas reservoirs. This procedure is similar to that 
followed by [Chapman et al.| ( |2004| l when measuring CO sizes in 
sub-millimeter galaxies. Individual objects vary somewhat in the 
surface brightness sensitivity achieved. Interpolating from the sur- 
face brightness profiles in Section |4] suggests that measuring the 
extent at this range of surface-densities introduces an error of upto 
~15% into the measured sizes. Given the small radial extent of the 
gas in many of our early-type objects, this error is likely insignif- 
icant when compared with the uncertainty due to beam smearing. 
The use of a scale radius obtained from a parametric fit to the CO 
surface density profile would be preferable, but the limited CO spa- 
tial extent (compared to the beam) in many objects, and the absence 
of an obvious suitable parametric form (e.g. Section |4| make this 
impractical at this point. 

The beam size of interferometric observations smears the 
edges of molecular gas distributions, making them appear larger 
than they truly are. As a first-order attempt to remove this beam 
smearing, the beam size in the direction of maximum extent was 
quadratically subtracted from the measured diameter, to give a 'de- 
convolved' diameter (d ) as follows: 

d = ^/d' - (2r-boam)2, (2) 

where 

^max^min 

^bcam , 1 

and 

db = 0GPA — Sbpa (4) 

Here d is the measured CO diameter, 6max is the beam semi-major 
axis, &min is the beam semi-minor axis, 6bpa is the beam position 
angle in degrees, (^gpa is the measured position angle along the 
maximal CO extent in degrees, rbcam is the radius of the beam at 
angle ^gpa, and 9], is the measured angle from the beam axis, to 
the molecular gas position angle. 

Figure [T] shows three examples of the redshifting and CO ex- 
tent measurement process. In some cases, such as NGC7331 shown 
in the top panel, the structure and extent of the molecular gas in the 
original map is preserved. In other cases, like IC342 shown in the 
middle panels of Figure [T] the extent of the molecular gas is simi- 
lar in the redshifted map, but the structure observed originally has 
disappeared. Finally, in some cases the spiral galaxies from BIMA- 
SONG would not be detected at all at the ATLAS^° distance and 
resolution. Such an example is NGC5247, shown in the bottom pan- 
els of Figure[T| 

In total, 3 1 of the 44 BIMA-SONG late-type galaxies had mea- 



surable CO extents after redshifting. The majority of the galaxies 
which are undetectable after redshifting become so because of the 
greater noise introduced. Galaxies with zero-spacing information 
available all have measurable extents after the redshifting proce- 
dure. Only the brightest 24 spiral-galaxies were selected to have 
on-the-fly maps created, and hence it is not clear if these galaxies 
survive redshifting better simply because they are brighter, or be- 
cause more of the large scale emission has been recovered (reduc- 
ing the effect of beam-dilution). Removing the four objects without 
zero-spacing information does not change our results, suggesting 
any bias introduced is small. 

40 of our 41 ETGs are included here. One ETC, NGC2685 
was removed from this analysis, as it has multiple components 
making measurement of the extent difficult. The field of view of 
the ATLAS''^ observations is often smaller than that achieved by 
BIMA-SONG. However, we do not expect this to bias our results, 
as very few of the ETGs have gas extended enough to reach the 
half-power beam point, and none reach to the edge of the field 
of view. The molecular gas extents of the ATLAS^"^ and BIMA- 
SONG galaxies, and the parameters used to derive them are listed 
m Tables[I]and|2] For the ATLAS'^° galaxies, the beam position an- 
gle and the beam major and minor axes used in the beam correction 
procedure are taken from the original data source, as listed in Table 
□ 

3.3 Comparing spirals and ETGs 

The distribution of absolute molecular gas extents in the redshifted 
BIMA-SONG and ATLAS^°galaxies are plotted in Figure H] The 
average BIMA-SONG galaxy CO radial extent is 2.1 kpc, while 
that for ATLAS^"^ galaxies is 1.1 kpc. The average surface density 
these are measured at is the same (Eh2 = 15 Mq pc~^). A Mann- 
Whitney U (MW-U) test gives only a 3% chance that they are drawn 
from the same parent distribution. We prefer the Mann- Whitney U 
test to a Kolmogorov-Smimov (KS) test here, as it is less sensitive 
to differences in the shapes of the distributions in the two samples, 
which are uncertain due to the low number of mapped objects. In 
all cases presented here, using a KS test gives similar results. 

We are also able to estimate how much more extended the 
ATLAS^° galaxies would have to be, on average, for the distribu- 
tions to be consistent. We do this by finding a multiplicative scaling 
factor that maximises the statistical similarity between the distribu- 
tions (as calculated from a MW-U test). Performing this analysis, 
we find a scaling factor of 1.6. The measured average surface den- 
sity increases by a factor of between 3 and 1.5 when one decreases 
by a factor of 1.6 the radius at which the CO extent is measured, 
suggesting that the molecular gas really is much more extended (in 
absolute terms) in late-type galaxies than in ETGs. 

The average distribution of stars in an ETG and in a spiral 
are quite different, the former generally being more concentrated 
because of the contribution of a large bulge. Indeed we find here 
that the effective radius of the BIMA-SONG galaxies is 2 kpc larger 
than the ATLAS'^'^ galaxies, and thus it is important to investigate 
the gas extent relative to the extent of the stellar matter. The top 
panel of Figure[3]shows the molecular gas extent normalized by the 
effective radius (i?c) of each galaxy. This is the standard scalelength 
for ellipticals, and represents the radius at which half the light is 
contained, assuming circular symmetry. The values used here 
come from Paper I for ATLAS^° galaxies, and for BIMA-SONG 
we use a combination of the RC3 ( |de Vaucouleurs et al.|199l) and 
2MASS i jjarrett et al.|2000| measurements. In an identical way 
as was done in the ATLAS'^'' sample paper, we use a scaled mean 
of the RC3 and 2MASS Rc values, where the scaling factor takes 
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Figure 1. Three example CO(l-O) moment zero maps from the BIMA- 
SONG survey, presented as observed in the left column, and redshifted us- 
ing the procedure outlined in Section [3?T] in the right column. The images in 
the left column are 85" square. The angular scale on the right hand images 
varies dependant on the distance to the object. The difference in beam size 
(displayed as a black oval in the lower left of each plot) gives an indication 
of the relative distance, and the new angular scale. The arrow in the upper 
and middle panels on the right indicate the measurement of maximum extent 
before beam subtraction. In the bottom right figure no emission is detected 
above 3o". 



into account the wavelength dependance of the effective radius for 
a given stellar population (see Paper I for the full justification). 



- 0.5 (i?o,H,C3 + ^2MASS X i?e, 



2MASSJ, 



(5) 



where i?e,RC3 is the RC3 value at B-band, i?c,2MASS is the me- 
dian of the J,H and Ks-hanA values in the 2MASS catalogue. 
The X2MASS factor depends on the average stellar population of 
the galaxies. Paper I showed a X2MASS of ~1.7 is appropriate for 
ETGs. Using the 2MASS and RC3 values for the BIMA-SONG 
spirals (in the same way as Paper I) we calculate that an X2MASS 
factor of ~1.26 is more appropriate in these late-type galaxies. We 
note however that if we used the same X2MASS factor for both sets 
of galaxies this would not change our results. 

The top panel of Figure[3]shows that the average BIMA-SONG 
galaxy normalized radial extent is 0.54, while the ATLAS'^'^ galaxy 
normalized mean is 0.47, however this difference is not statistically 
significant. A MW-U test shows we are unable to reject the hypoth- 
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Figure 2. Histogram showing the absolute molecular gas extent of the 
BIMA-SONG (blue) and ATLAS^^ (red) sample galaxies. The ATLAS^'^ 
histogram contains 40 galaxies, and the BIMA-SONG histogram contains 
31. The dashed lines show the mean extent of the galaxies in that histogram. 



esis that they are drawn from the same parent distribution (with a 
probability of 46%). A similar analysis to that conducted above sug- 
gests that the distributions are close to maximum similarity, with 
a scaling factor of 1.02. Overall this suggests that on average the 
molecular gas in early- and late-type galaxies covers a similar frac- 
tion of the effective radius. 

In an attempt to verify if the reduction in the average extent 
difference is due to the normalization radius adopted, we next at- 
tempted to normalize the CO diameter by i?25, the radius of the 
25 mag arcsecond"'^ isophote in B-band. We use the average R25 
value from LEDA ^Paturel et al. [19911 for each galaxy. The results 
are shown in the middle panel of Figure |3] The average BIMA- 
SONG galaxy normalized radial extent is 0.18, and the ATLAS^° 
normalized mean is 0.16. The small extent difference implied is not 
however statistically significant. A MW-U test gives a 28% chance 
that the BIMA-SONG and ATLAS^° normalized extents are drawn 
from the same parent distribution (a scaling factor of 1.14 is re- 
quired to maximise this probability). On average, once again, the 
molecular gas in these systems seems to cover a reasonably similar 
fraction of the 25 mag arcsecond"^ isophote. 

The result that the molecular gas distributions have similar dis- 
tributions of CO extent to the scale radii suggests that on average 
the molecular gas in larger galaxies is more extended, and that there 
is little dependence on Hubble type. If this is true, one would ex- 
pect a similar trend when normalizing by the i^^-band luminosity 
(bottom panel of Figure|3|. The Lk, values used for BIMA-SONG 
galaxies are the average value from the LEDA catalogue ( |Paturel| 



et al. 



1991 1 



and for ATLAS'''^ galaxies we use the value tabu- 
lated in Paper I. These absolute magnitudes are converted into a 
luminosity assuming that the absolute magnitude of the sun at Ks- 
band is 3.28 mag (Binney & Merrifield|199 8l. The average BIMA- 



SONG galaxy normalized radial extent is 10 
the ATLAS'^'^ galaxies normalized mean is 10" 



kpc Lq , while 
kpc Lg^. A 

MW-U test gives a 40% chance they are drawn from the same par- 
ent distribution, and this distribution is close to maximal similarity 
(scaling factor of 0.997). 

The significance of this similarity between the extent of the 
gas in ETGs and spirals in the context of their evolution will be 
discussed further in Section [6.1.1l 
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Figure 3. As Figure |2] but showing the relative molecular gas extent of the 
BIMA-SONG and ATLAS'^'^ sample galaxies. This has been normalized 
by the scaled mean of the effective radius of the galaxy al B, H, J and Kg 
band in the top panel, by the radius to the 25 mag arcsecond^^ isophote 
in B-band in the middle panel, and by the galaxies i^^-band luminosity in 
the bottom panel. The dashed lines show the mean normalized extent of the 
galaxies in that histogram. 



3.4 Effect of environment and gas origin 

We can also break down the ATLAS'"^ sample to investigate trends 
in the molecular gas extent with other parameters. The top panel 
of Figure [4] shows the extent of the molecular gas in the ATLAS'^'^ 
galaxies normalised by the effective radius, split into those galaxies 
where the kinematic misalignment between the molecular gas and 
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Figure 4. The relative molecular gas extent of ATLAS^'^ sample galaxies, 
normalized by the galaxies effective radius. In the top panel the ATLAS^'^ 
galaxies are split into those with kinematically-aligned molecular gas (top), 
and misaligned (bottom) with respect to the stars. Kinematic misalignments 
are taken from Paper X. In the bottom panel the ATLAS^'^ systems are split 
into Virgo cluster members (top) and field galaxies (bottom). Virgo cluster 
membership is taken from Paper I. 



the stars is greater than 30 degrees (see Paper X), and those with 
kinematically-aligned molecular gas. The aligned galaxies have an 
average extent of O.SORo and the misaligned galaxies 0.58 Rc - The 
small difference between these figures is not statistically significant, 
as a MW-U test gives a 26% chance that the kinematically-aligned 
and misaligned galaxies are drawn from the same parent population 
(a scaling factor of 1 .23 maximizes this probability). 

In Paper X, we also reported that the origin of the molecu- 
lar gas displayed a strong environmental dependence. The bottom 
panel of Figure [4] shows a comparison of the extent of the molecu- 
lar gas, normalised by the effective radius, in the Virgo cluster and 
field ETGs. The Virgo cluster galaxies have molecular gas distribu- 
tions that extend on average to 0.39 Rc, while field galaxies have 
more extended molecular gas distributions, reaching out to 0.60 Re. 
A MW-U test gives only a 3% probability that the Virgo and field 
ETGs are drawn from the same parent distribution, and a scale fac- 
tor of 1.87 is required to maximize this probability. We discuss pos- 
sible reasons for this environmental difference in Section l6.1.2l 

Given this environment dependance we revisited the correla- 
tions for spirals and ETGs discussed above, this time comparing 
only the field ETGs with the BIMA-SONG spirals (all of which are 
in the field). Removing the cluster galaxies leads us to estimate the 
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mean absolute extent of the gas in field ETGs as 1.23 Kpc. Thus 
the average field spiral galaxy has gas which is ~1.5 times as ex- 
tended as the average field ETG. When considering the normalised 
trends however, the increase in the average normalised radius actu- 
ally makes it statistically more likely that the field ETGs and spirals 
are drawn from the same population, strengthening our conclusions. 

4 Gas Surface Brightness Profiles 



Various authors (e.g. |Young & Scovillel|1982[ |Regan et al.| 
|2001[[Leroy et al.|2009[ l have studied the correspondence between 
the radial profile of the molecular gas and that of the stars in spiral 
galaxies. On small scales (~100 pc), the clumpy nature of the ISM 
causes the molecular gas to deviate from a simple exponential pro- 
file, but at larger scales (such as those probed by our CARMA ob- 
servations; 4'.'5 at 24.5 Mpc corresponds to ~ 500pc) the molecular 
gas in the disk of many spiral galaxies is found to have a reasonably 
similar profile to that of the stellar surface density, and has a similar 
scale-length. Although molecular rings, and excesses and deficits 
of central molecular emission are found in spiral galaxies, outside 
the bulge region the surface-brightness profile is still often compa- 
rable to the stellar surface brightness. As the overall distribution of 
extents is similar in spirals and ETGs, one might expect a similar 
result, with the molecular gas density profiles following those of the 
stars. However, the observed molecular gas structures in our galax- 
ies vary widely, and include central discs, rings, bars, spirals and 
disturbed gas. These galaxies are also much more bulge dominated 
than spirals. Molecular gas in ETGs could therefore have very dif- 
ferent surface brightness profiles than the relaxed discs seen in the 
majority of spiral galaxies. 

In order to investigate this issue, we extracted azimuthally av- 
eraged radial surface brightness profiles from the undipped CO in- 
tegrated intensity maps (to avoid any biases against faint emission) 
for all the sample galaxies. We used the CO kinematic position an- 
gle (from Paper X) and the 'best' inclination for the molecular gas 
(from Table 1, column 8 in Paper V) to estimate the average CO 
brightness in concentric elliptical annuli of one beam width (cen- 
tred at the optical nucleus of the galaxy as tabulated in Paper I). 
For the galaxies close to edge-on (inclination >80°), as we do not 
resolve the molecular gas discs thickness we estimate the surface 
brightness in rectangular regions one beam width thick along the 
major axis of the disc. These CO surface brightness profiles were 
then compared with the r-band stellar luminosity surface density 
profile convolved to the same resolution. The stellar surface bright- 
ness profiles was extracted from the MGE (Emsel lem, Monnet &| 



Bacon 



I994| model of each ATLAS^° galaxy (Scott et al., in prep) 



with the software of |Cappellari| ( |20()2) . We used the MGE values 
rather than those extracted directly from the r-band images, as the 
MGE models have been carefully fitted to remove the effects of 
foreground stars and other observational effects. MGE models were 
not available for NGC 4292 and PGC 0581 14 (due to a lack of suit- 
able r-band images), so in these two cases we used the method de- 
scribed above directly on A's-band images. Using Jfa-band images 
for all sources would not affect our results. 

In Figure |5] we show four examples of profiles chosen to be 
representative of the types of radial molecular gas surface density 
profiles we see. Profiles for all of the CO mapped ATLAS^° galax- 
ies are shown in Appendix A (in the online material). Some galaxies 
(such as PGC 0581 14, shown in the top-left panel of Figure |5} fol- 
low the stellar luminosity surface density closely, falling off with a 
similar scale length as the stars. We denote such systems as regu- 
lar in this work. We find that unlike in spirals where they make up 



~80% of the population jRegan et al.|200ll l, regular systems make 
up only half (21/41; 51%) of our ETG sample. 

Around a third (12/41; 29%) of the galaxies in this work have 
molecular gas radial profiles that fall off slower than those of the 
stars, and 7/12 of these systems appear to be truncated at large radii 
(the CO surface brightness does not smoothly drop to our detection 
limit, but ends abruptly). In each plot we include the first unde- 
tected radial bin with a 3a upper limit, to show this. The galaxy 
(NGC 3489) shown in the upper-right panel of Figure |5] is an ex- 
ample of such a profile. We classify such systems as having excess 
emission (with or without truncation) from now on. 

Approximately 7% of systems (3/41) have profiles that show a 
lack of emission in the central regions, with a peak further out be- 
fore the profile falls off (and is sometimes truncated). We have in- 
cluded the Helix galaxy (NGC 2685) in this category, based on the 
information available in |Schinnerer & Scoville|p002| l. The width of 
such features is normally about that of the synthesised beam, sug- 
gesting that they are not resolved. NGC 4324 (FigurejS] bottom left 
panel) is an example of such a profile (which we denote as rings). 

The last major category of profiles we see are like NGC 4710 
(bottom-right panel of Figure [5]l. The gas surface density falls 
quickly from a central maximum before peaking again at larger 
radii. Such systems make up ~7% (3/41) of the galaxies in our 
sample. These objects are all highly inclined, and such profiles are 
thus not present in the BIMA-SONG sample, which was selected 
to avoid edge on galaxies. We denote these systems as composites, 
and discuss their nature (and the nature of the other classes defined 
above) in Section|6j2] 

Two galaxies (NGC 1266 and UGC 05408) were not resolved 
at the resolution of our D-array CARMA data (but see |Alatalo et al.| 
|20Il 'for a full discussion of the gas in the molecular outflow of 
NGC 1266). We simply denote these surface-brightness profile as 
unresolved. Table |3] lists the surface brightness profile class as- 
signed to each galaxy. For a discussion of how these profile types 
arise, and any dependance on galaxy properties, see Section [6^ 

5 Gas liinematics 

5.1 Comparing CO velocities to circular velocities 

In Paper V of this series, we compared the CO single-dish and in- 
terferometric line-widths with JAM models of the circular velocity 
of each galaxy. We showed that most of the galaxies with double- 
peaked or boxy spectra are likely to have relaxed gas reaching be- 
yond the peak of the rotation curve. This allowed us to use these sys- 
tems to investigate the CO Tully-Fisher relation of ETGs. It is how- 
ever possible to go beyond this analysis, and for each galaxy com- 
pare the position-velocity diagram (PVD) of the mapped molecular 
gas with the circular velocity curve derived from the JAM models, 
as well as the observed stellar and ionised gas velocities (as done 
for four of these galaxies in |Young, Bureau & Cappellari|2008| . 

We extracted the CO PVDs from our data cubes along the 
kinematic position angle determined in Paper X, as described in 
A12. Two example PVDs are shown in Figure|6] while those for the 
whole of the sample are shown in Appendix B (in the online mate- 
rial). We do not include the PVDs for the two non-ATLAS^° galax- 
ies (no corresponding IFU data), NGC 2685 (kinematics discussed 
in |Schinnerer & Scoville|2002^ , NGC 1266 (molecular outflow dis- 
cussed in |Alataloetal.|201I[ >, PGC0581 14 (insufficient photometry 
for a JAM model) and UGC05408 (which is unresolved at the reso- 
lution of our data). 

The observed PVDs are overlaid with the projected JAM cir- 
cular velocity curves in black (see Section [23] for details). If the 
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Table 1. Parameters used in calculating the maximal extent of the molecular gas in the ATLAS'^'^ early-type sample galaxies. 



Galaxy 


R25 


Deo 


^'gpa 


Rco 


«co 




log(^) 


Ref. 




(kpc) 


(") 


(deg) 


(kpc) 




(kpc/LQ) 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


1C0676 


8.81 


22.0 


164 


1.22 


0.46 


0.14 


-10.13 




1C0719 


5.39 


26.8 


50 


1.85 


0.48 


0.34 


-10.13 




1C1024 


3.27 


25.0 


23 


1.21 


1.08 


0.37 


-9.84 




NGC0524 


9.36 


19.4 


47 


1.05 


0.21 


0.11 


-11.17 




NGC1222 


5.18 


25.0 


72 


1.90 


1.00 


0.37 


-10.09 




NGC1266 


5.92 


8.9 


145 


0.62 


0.22 


0.10 


-10.65 




NGC2697 


8.07 


30.0 


58 


2.02 


0.64 


0.25 


-10.08 




NGC2764 


8.94 


36.5 


21 


3.42 


1.45 


0.38 


-10.05 




NGC2768 


23.72 


13.0 


20 


0.64 


0.10 


0.03 


-11.38 




NGC2824 


4.61 


16.0 


142 


1.36 


1.03 


0.29 


-10.31 




NGC3032 


4.40 


14.0 


315 


0.50 


0.54 


0.11 


-10.41 




NGC3182 


9.52 


16.5 


120 


1.16 


0.34 


0.12 


-10.52 




NGC3489 


6.05 


20.0 


44 


0.54 


0.43 


0.09 


-10.77 




NGC3607 


14.84 


30.0 


122 


1.49 


0.36 


0.10 


-11.03 




NGC3619 


10.47 


11.1 


85 


0.46 


0.13 


0.04 


-11.08 




NGC3626 


7.65 


27.0 


148 


1.22 


0.50 


0.16 


-10.55 




NGC3665 


11.85 


21.0 


65 


1.58 


0.32 


0.13 


-11.08 




NGC4119 


4.90 


17.6 


107 


0.63 


0.20 


0.13 


-10.55 




NGC4150 


3.98 


30.0 


161 


0.91 


0.78 


0.23 


-10.01 




NGC4292 


3.72 


9.2 


44 


0.31 


0.15 


0.08 


-10.44 




NGC4324 


6.63 


48.0 


48 


1.90 


1.19 


0.29 


-10.08 




NGC4429 


13.53 


14.0 


90 


0.48 


0.14 


0.04 


-11.36 




NGC4435 


6.71 


10.5 


23 


0.28 


0.11 


0.04 


-11.39 




NGC4459 


9.59 


20.0 


113 


0.51 


0.19 


0.05 


-11.16 


2 


NGC4476 


4.41 


21.6 


24 


0.78 


0.60 


0.18 


-10.13 


1 


NGC4477 


9.15 


7.4 


47 


0.21 


0.07 


0.02 


-11.50 


5 


NGC4526 


17.52 


14.0 


111 


0.44 


0.13 


0.03 


-11.51 


2 


NGC4550 


7.48 


8.0 


1 


0.30 


0.26 


0.04 


-10.74 


4 


NGC4694 


7.61 


26.0 


97 


1.01 


0.43 


0.13 


-10.17 




NGC4710 


6.25 


60.0 


30 


2.41 


0.96 


0.39 


-10.34 




NGC4753 


19.13 


29.9 


92 


1.59 


0.29 


0.08 


-11.15 




NGC5379 


9.57 


41.3 


57 


2.95 


0.97 


0.31 


-9.67 




NGC5866 


10.16 


63.0 


144 


2.26 


0.86 


0.22 


-10.56 




NGC6014 


9.18 


9.7 


128 


0.46 


0.12 


0.05 


-10.83 




NGC7465 


4.36 


25.6 


71 


1.67 


1.75 


0.38 


-10.17 




PGC029321 


4.26 


11.3 


125 


0.84 


0.55 


0.20 


-10.05 




PGC058114 


2.07 


16.1 


91 


0.77 


0.72 


0.37 


-10.05 




UGC05408 


6.80 


5.9 


120 


0.66 


0.43 


0.10 


-10.31 




UGC06176 


7.88 


8.0 


36 


0.32 


0.16 


0.04 


-10.87 




UGC09519 


3.10 


12.9 


145 


0.69 


0.77 


0.22 


-10.21 





Notes: Column 1 lists the names of the galaxies in the ATLAS'^'-' sample considered here. Column 2 lists the radius of the 25 mag arcsec"^ 
isophote at B-band, extracted from HyperLEDA and converted into kpc using the distance to the galaxy as tabulated in Paper 1. Column 3 
contains the maximum CO angular diameter, calculated as described in Section [33] Column 4 contains the galaxy position angle. Column 5 
contains the maximum CO extent, beam corrected and converted to a linear size using the distance to the galaxy. Column 6 is the ratio of the 
CO extent to the effective radius of the galaxy (tabulated in Paper I). Column 7 lists the ratio with respect to R25 . Column 8 lists the ratio to the 
Ks-band luminosity of the galaxy, calculated from the magnitudes listed in Paper I, and converted into a luminosity assuming that the absolute 
magnitude of the Sun at Ks-ba.nd is 3.28 mag jBinney & Merrifield 19981. Column 9 lists the data reference for the CO interferometry. The 
beam position angle and the beam major and minor axes used in the beam correction are taken from these papers: (1) A12 (2) [Young, Bureau] 
[& Cappellari]j2008f , (3) [Crocker etarH20O8f , (4) [Crocker etarH2009f , (5) [Crocker eFarHMot . 



molecular gas is counter-rotating, we also overlay a mirrored ver- 
sion of the circular velocity curve, to allow a direct comparison with 
the CO, and to investigate if they are in agreement. If the CO tightly 
follows the circular velocity curve, this suggests that CO is dynam- 
ically cold and relaxed, even in ETGs. Substantial disagreement 
would indicate that the gas is dynamically disturbed, not on cir- 
cular orbits, or indicate problems with the JAM models (e.g. invalid 
assumption of a constant anisotropy or mass-to-light ratio (M/L); 
see |Young, Bureau & Cappellari 2008| l. 

By inspection of the PVDs and circular velocity curves, we es- 
tablish that the molecular gas reaches beyond the turnover radius in 



68±6% of our ETGs (as indicated in Table|3j. This compares well 
to the 71 ±7% of galaxies which were classified as having double- 
peaked or boxy single-dish spectra in Paper V (and hence were as- 
sumed to reach beyond this turnover). 

5.2 Simulated observations 

To investigate if the observed PVDs are consistent with dynam- 
ically cold gas rotating at the predicted circular velocities, while 
including observational effects, we have developed a new IDL ob- 
servation simulation tool. This routine, the KINematic Molecular 
Simulation (KinMS) tool, allows one to create a mock interfero- 
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Table 2. Parameters used in calculating the extent of the molecular gas in the BIMA-SONG late-type galaxies. 



Galaxy 


Dist. 




Rc 


i?25 


Deo 


0GPA 


Sbpa 






Rco 


«co 

-Ro 


«co 
-R25 






(Mpc) 


(mag) 


(kpc) 


(kpc) 


(") 


(deg) 


(deg) 


(") 


(") 


(kpc) 






(kpc/L0) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


(14) 


IC342 


3.9 


-23.40 


9.80 


9.14 


10.3 


10 


33 


18.9 


17.2 


0.51 


0.05 


0.06 


-10.96 


NGC0628 


7.3 


-23.14 


4.42 


11.13 


34.3 


338 


8 


7.2 


5.3 


2.02 


0.46 


0.18 


-10.26 


NGC1068 


14.4 


-24.96 


3.10 


14.85 


49.7 


12 


16 


8.9 


5.6 


2.91 


0.94 


0.20 


-10.83 


NGC2903 


6.3 


-22.96 


3.11 


7.11 


39.9 


30 


21 


36.5 


34.1 


2.12 


0.68 


0.30 


-10.17 


NGC3184 


3.6 


-20.57 


2.32 


3.74 


3.3 


10 


11 


16.7 


15.2 


0.20 


0.08 


0.05 


-10.25 


NGC3351 


10.1 


-23.36 


3.68 


7.89 


11.5 


5 


6 


17.9 


12.6 


0.23 


0.06 


0.03 


-11.29 


NGC3368 


10.1 


-23.70 


3.03 


9.93 
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Notes: Column 1 lists the names of the galaxies in the BIMA-SONG survey that had measurable extents after undergoing the redshifting 
procedure outlined in Section|3J] Column 2 contains the distance to the galaxy. Column 3 is the absolute magnitude of the galaxy at Kg-hand, 
extracted from HyperLEDA. Columns 4 and 5 list the effective radius and the radius of the 25 mag arcsec"^ isophote at S-band, extracted 
from HyperLEDA and converted to kpc using the (real) distance to the galaxy. Column 6 contains the maximum CO angular diameter (at the 
new redshifted distance), calculated as described in Section [T2] Columns 7 to 10 contain the galaxy position angle, the beam position angle 
and transformed beam major and minor axes, respectively. Column 1 1 contains the maximum CO extent, beam corrected, and converted to a 
linear size using the new distance to the galaxy. Column 12 is the ratio of the CO extent to the effective radius of the galaxy (tabulated in Paper 
I). Column 13 lists the ratio of the CO extent with respect to -R25. Column 14 contains the ratio of the CO extent with respect to the Xs-band 
luminosity of the galaxy, calculated from the absolute magnitudes, which have been converted into a luminosity assuming that the absolute 
magnitude of the sun at Kg-banA is 3.28 mag ( |Binney & Merrifield|1998j l. 



metric data cube based on arbitrary surface-brightness and velocity 
profiles, with a realistic treatment of projection, disk-thickness and 
gas velocity dispersion, as well as simulating observational effects 
such as beam-smearing, velocity binning, etc. We make this tool 
freely available to the community, along with full documentation 
and examples, at |http://www.eso.org/^tdavis| 

In this work we set up the KinMS tool to simulate our C ARMA 
observations, imposing a fixed velocity resolution (10 km s^^ chan- 
nels), bandwidth (~900 km s^^), pixel size (1") and the angular 
resolution of our CARMA data (typically 4'.'5, A12). We make the 
assumption that the gas is dynamically cold, with a velocity dis- 
persion (Tgas ~ 8 km s~^, and that it rotates at the predicted JAM 
circular velocity. We give the molecular gas disk/rings a fixed ver- 
tical scale height of lOOpc, similar to that observed in the Milky 
Way (e.g. |Nakanishi & Sofue|2006l using any sensible value of this 
parameter does not change our results). We project the model to the 
observed inclination, as derived from the JAM models (see Section 



|2.3^ . We then extract a PVD in the same way as for the real data, as 
described in Section [5T1 

Given the above setup, each model is only left with two free 
parameters, controlling the scaling and shape of the surface bright- 
ness profile. For the majority of galaxies, we assume that the CO 
is distributed in an exponential disc, with a scale length varied 
to match our observations. In the systems with molecular rings, 
we characterised the surface brightness profile with a Gaussian of 
FWHM ~100 pc centred on the observed ring radius. The parame- 
ters and profile type for each galaxy are listed in Table [3] It would 
be possible to remove this modelling step by using the observed 
surface-brightness profiles presented in Section |4] or even the ob- 
served moment zero maps. However these observed profiles would 
have to be de-projected, and de-convolved in some manner to re- 
move the effect of the beam, both processes which introduce sig- 
nificant uncertainties, especially given the limited CO spatial extent 
(compared to the beam) in many objects. In this paper we only aim 
to assess the agreement between the JAM circular velocities and 
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Figure 5. Elliptically averaged radial surface brightness profiles of the stars (red dashed line; r-band) and molecular gas (blue solid line with errors) for four of 
our ATLAS'^'^ ETGs (normalized at the first CO datapoint). The galaxies shown are PGC058114 (top-left; classed as a regular' galaxy), NGC 3489 (top-right; 
classed as an excess galaxy), NGC 4324 (bottom-left, classed as a ring) and NGC 4710 (bottom-right; classed as a composite system). The CO beamsize is 
shown as a black dotted gaussian. The error bars on the CO measurements denote the width of the elliptical annuli in the x-direction and the RMS noise in 
the elliptical annulus in the y-direction. The letter(s) after the galaxy name denotes the profile class. "Reg" denotes a regular' profile, "E" an extended profile, 
"E(T)" an extended profile with truncation, "R" a ring, "C" a composite, and "U" an unresolved profile. 



our data, not the internal details of gas distribution, and using sim- 
plified surface brightness profiles suffices for this purpose. Tests on 
well resolved galaxies from this sample showed that using the ob- 
served profiles would not change our conclusions. We do however 
note any major discrepancies between the model and the observed 
PVD in Section [63] Two example models are shown in the right 
column of Figure [6] and models for the whole sample are shown in 
Appendix B (in the online material). 

Comparing the CO PVDs to the JAM model circular velocities 
and modelled PVDs (by eye) allows us to identify different classes 
of galaxies. These are described in detail below: 

• "G" (good) - In 51% (18/35) of our sample the models match 
well both the shape and position of the CO PVDs. 

• "O" (over) - In a further 20% (7/35) of the galaxies the shape 
of the JAM circular velocity curve and modeled PVD is similar to 
that of the observed CO PVD, but the models predict higher veloc- 
ities than those seen in CO. 

• "X" (disturbed) - In 1 1% (4/35) of the galaxies where the CO 
is disturbed (as identified in A 12), the CO appears to have little 
coherent rotation and its kinematics are, as expected, unrelated to 
the JAM model predictions. 

• "P" (polar) - In 9% (3/35) of the galaxies the molecular gas 



rotates around the polar plane of the galaxy (Paper X). The model 
PVDs are systematically above the observed rotation, as expected 
if the potential of the galaxy is flattened. 

• "?" (misc) - Of the remaining 9% (3/35) of galaxies, two 
(NGC 3032, NGC 2824) show signs that the assumption of a con- 
stant M/L as a function of radius is invalid (likely due to strong 
star-formation in the central parts), and one (NGC 0524) is so face- 
on that a large degeneracy exists between the model inclination and 



M/L (as discussed in Cappellari et al. 2006 1 



Membership of these classes was determined by comparing 
the model and observed PVDs, and noting significant discrepancies 
by eye. It is thus important to note that these classifications rely 
on the projected velocity differences between the model and ob- 
served PVD. Where the gas PVD is X shaped we only compare the 
model to the arm of the X-shape with larger extent, as in the inner 
arm is likely to come from gas in the X2 orbits within bars, which 
can be highly elliptical in shape (e.g. [Bureau & Athanassoula|1999[ 
[Athanassoula & Bureau] 1999) . These galaxies almost all also show 
evidence for bars in the stellar light profile/kinematics (see Section 

In order to more quantitatively compare the velocity of the CO 
at a given radius with other measures of rotation we extracted a 
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Figure 6. CO Position- Velocity diagrams (orange with black contours) of NGC 3607 (top) and NGC 4429 (bottom), overlaid with the JAM circular velocity 
curve (thick black line). On the left, these have also been overlaid with the observed stellar (red dashed line) and ionised gas (blue line) rotation curve. On 
the right we overlay the modeled observations of each galaxy (blue contours), created assuming that the gas is dynamically cold and rotating at the predicted 
circular velocity, as described in Section [5^ These systems have good JAM models predictions of the circular velocity, as denoted by the profile class "G" 
displayed after the galaxy name. 



'trace' from the observed and model PVDs. This trace was con- 
structed by finding the velocity of the pixel with the highest CO flux 
at each radius in the PVD (as long as this pixel has a flux greater 
than three times the RMS noise in the cube). This trace will likely 
lie below the true circular velocity of the gas (as discussed in Sec- 
tion [63T|(, so some offset between this measure and the other ve- 
locity tracers is expected. Given this caveat however, it is possible 
to use this trace as a measure of relative variations between sample 
members. Edge on galaxies, where the line-of-sight passes through 
material at many different radii, could be expected to show system- 
atic offsets from less edge on systems. This effect is not observed 
in this sample however (the edge on systems do not lie show any 
systematic offset from the other systems in any analysis we have 
conducted), likely due to the proper treatment of the 3D structure 
of these systems within the KinMS code. 

In Figure |7] we plot the average equivalent width (EW) of the 
H/3 emission against the normalized average difference between the 
CO trace and the model trace. This measure allows us to estimate 
the fractional difference between the observed CO velocity and the 
model. The average EW is the H/3 luminosity weighted mean cal- 
culated of all SAURON bins within an elliptical aperture of the same 
extent, axis ratio and position angle (Table[T|l as the molecular gas. 
Black points identify galaxies where the JAM model prediction of 



the circular velocity is good (black stars indicate galaxies where the 
gas is counter-rotating; see Paper X). Red points identify galaxies 
in which the JAM model prediction is over the CO velocities. The 
red points, where the JAM model predicts higher velocities than 
observed in CO, although not selected via this diagram, have as 
expected larger offsets than the black points, with (de-projected) 
offsets of <-0.2 dex (equivalent to a higher prediction of the veloc- 
ity by ^30 km s^^). Importantly, however, almost all the galax- 
ies with EW(H/?) A have significant (<-0.3 dex) offsets, and 
most are classified as having JAM models that predict higher veloc- 
ities than those observed in CO. We discuss these findings, and why 
these classes arise in Section l6.3.1l 

It is clear that our discrete classification by eye of the projected 
PVDs is not completely clean, and a few galaxies could be put in ei- 
ther the good or over category. The exact fraction of galaxies in the 
good and over classes are uncertain at the ~5-10% level, however 
this is not enough to change the results presented here. 

5.3 Stellar velocities 

In the left column of Figure[6]and Appendix B (in the online mate- 
rial) we also overplot the observed stellar rotation curves onto the 
CO PVDs with red dashed lines. These were extracted from the 
SAURON IFU data, using pseudo-slits along the same PAs as those 
used for the CO (tabulated in Paper X). The observed velocity field 
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Figure 7. Normalized radially averaged velocity differences between the 
trace of the CO and model PVD, plotted against the luminosity average 
equivalent width of the H/3 emission within an elliptical aperture (with the 
same extent, PA and axis ratio as the molecular gas). Black points denote 
galaxies in which the JAM model prediction of the circular velocity matches 
well the observed molecular gas rotation (black and red stars indicate galax- 
ies where the gas is counter-rotating; see Paper X). Red points are those 
galaxies in which the JAM model predicts higher velocities than those ob- 
served in CO. All other JAM profile classes are not plotted on this diagram. 
The dashed line indicates a velocity percentage offset of zero. NGC5866 is 
missing from this figure because it lacks a measurement of the equivalent 
width of H/3. 
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Figure 8. Radially averaged stellar velocity dispersions predicted using the 
simple approximations in Equation 9, plotted against the observed stellar 
velocity dispersion (Paper I) within one effective radius. The error bars on 
(Tfle are not the formal error, but show the variation in observed velocity 
dispersion along the slit. The points in this figure are colour coded to denote 
the JAM model profile class. Red and black points have the same meaning as 
in Figure|7] and in addition red stars demote galaxies with polar molecular 
gas, open red squares represent disturbed galaxies and open blue squares 
are the remaining systems. Stars of these colours indicate galaxies which, 
in addition to the above criteria, also have misaligned molecular gas. The 
dashed line shows the one-to-one relation as a guide to the eye. 



was fitted with a three dimensional (minimum curvature sphne) sur- 
face, and the required values along the slit were interpolated. This 
process allows us to minimize the effect of bad bins which happen 
to lie within the slit. 

Stellar systems with large velocity dispersions rotate slower, at 
a given radius, than dynamically cold components like the molec- 
ular gas. The difference between the velocities is because the stars 
do not move on circular orbits (known as asymmetric drift in a disk 
system). One would thus expect the velocity difference between the 
stars and the molecular gas to be larger in higher dispersion (and 
thus likely more massive) galaxies. Comparing the stellar rotation 
to a cold gas rotation tracer in principle allows one to estimate the 
magnitude of the asymmetric drift in each galaxy and disentangle 
these from the effect of turbulence, using e.g. the formalism of |Wei-| 
[jmans et^ ( |2008| l. Such a detailed comparison is outside the scope 
of this work, but will be attempted in a future work in this series. 
Here we derive a simple zeroth order estimate of the velocity dis- 
persion in these galaxies by assuming that the second moments of 
the molecular gas and stars will be equal, as they trace the same 
underlying potential: 

VV.^ + ol = ^/Vi,, + ai^, (6) 
and thus 

Cr*, prod = \/(VgL-K2)+0-|as, (7) 

where V„ and Vg^s are the measured stellar and molecular gas ve- 
locities at a given point, o-gas is the molecular gas velocity disper- 
sion (here we assume 8 km s^^, the dispersion of HI in the solar 
neighborhood; [van der Kruit & Sho stak 1982; Dickey, Hans on &| 
|Helou|1990| >, and cr,^ prod is thus the predicted stellar velocity dis- 
persion required to produce the observed velocity differences be- 
tween the CO and stars. 



Figure [8] shows a comparison of the radially averaged pre- 
dicted velocity dispersions derived using the zeroth order approx- 
imation described above, plotted against the observed stellar ve- 
locity dispersion of each galaxy (Paper I) measured within one ef- 
fective radius. The error bars on aue are not the formal error, but 
show the variation in observed velocity dispersion along the slit. 
The points in this figure are colour coded to denote the JAM model 
profile class; red and black points have the same meaning as in Fig- 
ure|7] and in addition red stars denote galaxies with polar molecular 
gas, open red squares represent disturbed galaxies and open blue 
squares are the remaining systems (discussed in Section [5T| . The 
dashed line shows the one-to-one relation as a guide to the eye. 

Figure[8]shows that the velocity difference between the molec- 
ular gas and stellar velocities does appear to scale with the degree 
of pressure support in the system, as expected. The highest velocity 
dispersion systems (a > 200km s~^) show some sign of preferen- 
tially lying above the relation, which may indicate that the zeroth 
order approximations we use to derive the predicted velocity dis- 
persions break down as the velocity dispersion becomes similar to 
the circular velocity. Given the low number of objects, however, 
the statistical significance of any such trend is low. Figure [S] also 
highlights that the systems which show differences between the CO 
rotation and the JAM prediction (red solid points) tend to have low 
(>120km s~^) velocity dispersions. 

5.4 Ionized gas velocities 

In the left column of Figure |6] and Appendix B (in the online ma- 
terial), we also overplot the observed ionised gas rotation curves 
onto the CO PVDs with a blue line. These were extracted from the 
SAURON IFU data, using a pseudo-slit along the same PA as used 
for the stars and CO (tabulated in Paper X). The observed velocity 
field was fitted with a surface, and the required values along the slit 
interpolated, as was done for the stars. 

Ionized gas discs would normally be expected to rotate slower 
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Figure 9. Radially averaged velocity difference between the molecular and 
ionised gas plotted against the luminosity weighted EW(H/3) (top panel; 
as in Figure[7j and the [0III]/H/3 ratio (bottom panel). Symbols, and their 
colour coding have the same meaning as in Figure |7] The dashed line de- 
notes zero velocity offset 



than the dynamically cold molecular gas, due to their larger dis- 
persion, but faster than the dispersion-dominated stars. The relative 
velocities of the ionized gas and the CO can provide clues on the 
source and ionization mechanism of the gas emission. For instance 
strong H/3 emission in these molecular gas rich galaxies is often 
associated with HII regions embedded in the dynamically cold star- 
forming disk, while regions where [OIII] emission dominates seems 
to have a significant contribution from a dynamically hotter compo- 
nent of the ionized gas which is not related to star formation (e.g. 
ISarzi et al.|2"006l|Young, Bureau & Cappellari|2008l >. 

The top panel of Figure|9]shows the radially averaged velocity 
difference between the molecular and ionised gas plotted against 
the average EW(H/3) (as in Figure |7]l. The bottom panel shows the 
same measure plotted against the mean [OIII]/H/3 ratio, measured 
over the same area as used for the EW(H/5) measurements. Galax- 
ies with an EW(H^) > 3A and those with log([OIII]/H/3) < 0.0 
have ionised gas velocities close to those of the molecular gas. Sys- 
tems with a lower EW(H/3), and higher [OIII]/H/? ratios have a wide 
spread of velocity differences, and many have CO profiles which 
agree well with the JAM model prediction. These findings are dis- 
cussed in more detail in Section l6.3.3l 



Table 3. Surface brightness profiles 


and position-velocity diagram classes 


and model parameters for the ATLAS^'' ETG sample. 
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Notes: Column 1 contains the galaxy name. Column 2 lists the surface bright- 
ness profile class described in Section |4] "Reg" denotes a regular profile, "E" 
an excess profile, "E(T)" an excess profile with truncation. "R" a ring and "C" 
a composite profile. Column 3 lists the position-velocity class described in Sec- 
tion |5.2| "G" denotes a good profile, "O" a profile where the JAM model predicts 
a higher rotation velocity than observed in CO, "P" denotes a polar molecular 
gas disk, "X" denotes a galaxy with disturbed gas, and "?" denotes systems 
with known problems. A dashes shows that the galaxy was not included in this 
analysis, for reasons discussed in Section|5J] Column 4 shows the scale radius 
used to create the model PVDs described in Section Isl2l Column 5 shows the 
type of surface brightness profile used to create the model PVD: "e" denotes 
an exponential disc and "r" a Gaussian ring. Column 6 contains an "x" if the 
molecular gas extends beyond the turnover of the JAM circular velocity curve, 
and a dash if the gas is insufficiently extended (in two cases we do not have the 
data to determine this, these galaxies are denoted with a "o" in this column). 
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6 Discussion 

6.1 Molecular gas extent 

6.1.1 ETGs and spirals 

Earlier in this paper we identified scenarios that would result in 
the gas in ETGs being less extended than that in spiral galaxies. 
In Section |3.3| (Figure [2](, we indeed found convincing evidence 
that the absolute extent of the molecular gas in ETGs is smaller 
than that in spiral galaxies. However, spiral galaxies are generally 
more extended (due to their exponential light profiles; the average 
effective radius of the BIMA-SONG galaxies is 2 kpc larger than 
that for ATLAS^° galaxies). Pi gurejsjshows that the distribution of 
molecular gas extents seem to scale with the size and luminosity of 
the galaxy. That is, the distributions of the relative CO-extent (with 
respect to i?c, R25 and A's-band luminosity) of ETGs and spirals 
are similar. 

A similar correlation between the size scale of spiral discs and 
their large-scale CO distribution has been discussed before by sev- 
eral authors. |Young & Scoville] ( |1982| > found that the major- and 
minor-axis profiles of the molecular gas in NGC 6946 and IC 342 
were in good agreement with optical B-band radial profiles. [Reganj 
|et al.|p001^ found a similar result for a subsample of 15 BIMA- 
SONG spirals, for which the mean ratio of the CO to stellar (ex- 
ponential) scale lengths was 0.88=h0.14. [Leroy et al.|p009^ also 
find an average CO scale-length of 0.2±0.05 R25 for their HERA- 
CLES sample spiral galaxies. Although these parametric scale radii 
are hard to compare with the (isophotal) radii used in this paper, it 
appears that some link exists between the average size of the molec- 
ular reservoirs in galaxies and their stellar distributions. 

Despite the large intrinsic variation in scale-length between 
galaxies, [Regan et al.| ( |20 01 ) suggest that in a time-averaged sense, 
the stellar and molecular discs must be closely coupled. They posit 
that this coupling could arise from star formation associated with 
the molecular gas, that would over time cause the scale lengths to 
become equal. Alternatively they suggest that gravitational effects 
could exert some feedback on the molecular disc, affecting angular 
momentum transfer, damping out variations and causing the molec- 
ular disc to conform more closely to the stellar potential. Here we 
suggest that mass loss from previous generations of stars could also 
play a role. The mass loss material would interact, exchange an- 
gular momentum, and provide a means of coupling the stellar and 
molecular discs. 

The fact that we find a similar relative molecular gas extent 
in both the ATLAS^° ETGs and BIMA-SONG spirals leads to the 
suggestion that some combination of these processes that |RegMi| 
|et al.| ( [200T| l suggested were active in spirals could also be at work 
in ETGs. When one considers the low molecular gas mass fractions 
and star-formation rates present in our ETGs (typically molecular 
to stellar mass ratios for our ETGs are ~0.5%, as compared to 
>10% for spirals; e.g. [McGaugh & de Blok||1997^ , gravitational 
processes or stellar mass loss coupling seem more likely to be the 
dominant driver(s). Both of these processes require the gas to have 
been present and relaxed in the potential of the galaxy for some 
time, however. While the majority of the ETGs in our sample have 
dynamically relaxed gas (See Section|5j, often the gas is misaligned 
from the stellar body, and/or not distributed in a 'disky' way (Sec- 
tion |4|. As we find a similar extent for most subsamples of our 
ETGs (e.g. misaligned and aligned gas reservoirs; Section[3^, this 
suggests that the dominant process which controls the extent of the 
gas has to work on short timescales, and must be independent of the 
orientation of the gas within the galaxy potential. 



6.1.2 Gas kinematic misalignments and extent 

In Section [34] we investigated the extent of the molecular gas in 
subsamples of the ATLAS''^ ETGs. We found in the top panel of 
Figure|4]that kinematically-misaligned galaxies (where the gas has 
an external origin) have a similar distribution of extents as aligned 
galaxies (where the gas could be from internal sources). This re- 
sult is initially hard to understand, as gas with an external origin 
can have any amount of angular momentum, aligned in any direc- 
tion, and thus the radii at which incoming gas will settle depends 
strongly on the initial parameters of the source (e.g. [Barnes & HenT] 
|quist|1996l [Bois et al.|2010| l. In Paper X, however, we discussed 
evidence that the majority of the gas in (field) ETGs could be of 
external origin, even when it is aligned with the stars. If all the gas 
was from external sources, regardless of misalignment, this could 
provide a natural explanation for the similarity between the two 
distributions. Indeed, when one removes cluster galaxies (which 
we showed are more likely to have internally generated gas) the 
kinematically-aligned and misaligned field galaxies extent distribu- 
tions are still statistically indistinguishable (A MW-U test gives a 
30% chance that the BIMA-SONG and field galaxy ATLAS^° nor- 
malized extents are drawn from the same parent distribution). 

Molecular gas (and hence likely the star formation) in spi- 
rals, however, cover a similar proportion of the galaxy as discussed 
above. In a fast-rotator with a spiral progenitor, gas created from 
stellar mass loss (which comes predominantly from stars <1 Gyr 
old) would thus occupy a similar area. 

6.1.3 Gas extent and environment 

In Section [3^ we showed that the molecular gas in Virgo cluster 
ETGs covers only around half as much of the effective radius as 
the gas in an average field galaxy. This result is statistically sig- 
nificant at the 97% level, even given the low number of mapped 
galaxies in the Virgo cluster. This suggests that entering the cluster 
has significantly affected these molecular reservoirs. HI in cluster 
environments is known to be substantially affected by ram-pressure 
stripping, but it is thought that denser molecular gas can remain 
bound until much deeper in the inter-cluster medium (e.g. |Vollmer| 
|et al.|[2005| l. It is however possible that ram pressure could have 
compressed the molecular gas reservoirs in these objects. 

Paper X reports that all but one of these fast-rotating Virgo 
cluster galaxies have kinematically-aligned gas, and thus the molec- 
ular reservoir has either been substantially affected by the cluster 
environment or new gas has been generated internally. The mis- 
alignment of the gas in field ETGs was not found to be radius de- 
pendent at the resolution limit of Paper X. Removal of the outer 
gas, as implied here, would thus not naturally lead to the gas being 
measured as more aligned. 

Over time it is expected that stellar mass loss from old stars 
and any kinematically-misaligned molecular gas would collide, re- 
ducing the velocity and angular momentum of the gas, and causing 
the molecular reservoir to become smaller. If such a process is im- 
portant one would expect misaligned galaxies to show a negative 
correlation between the Ar parameter (a proxy for the amount of 
rotation present within the stellar component of the galaxy; see Pa- 
per III) and the extent of the misaligned gas reservoir Indeed, in 
our sample we do see a weak correlation in this sense, with faster 
rotating ETGs having less extended misaligned gas disks, but as 
we only have a total of twelve misaligned objects further work will 
be required to determine if this trend is confirmed in a statistically 
significant sample. 

Given the above mechanism, the observation that the 
kinematically-aligned gas in Virgo galaxies is less extended is con- 
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sistent with a scenario where some proportion of cluster ETGs at 
some stage had misahgned gas, which has been forced to align (due 
to the stellar potential and interaction with stellar mass loss). Inter- 
actions between misaligned molecular material and the ICM could 
also have a similar effect, although to first order one would expect 
the pressure exerted by the ICM to not correlate with the stellar 
rotation axis, and thus to create both co- and counter-rotating gas 
reservoirs. Any of these processes would have to force the gas into 
an aligned state quickly, however, as: 

(i) we observe very few systems in the cluster that still have 
kinematically-misaligned gas. 

(ii) the molecular gas fractions observed inside and outside the 
cluster are similar. Assuming the galaxies entered the cluster look- 
ing like field ETGs, and that they form stars with the same effi- 
ciency as field ETGs, then star-formation would use up significant 
fractions of the molecular material if the timescale for the gas to 
relax into the equatorial plane were long. 

If the molecular gas can be held stable against star formation (or 
equivalently have its star formation efficiency reduced) then the 
problem of gas depletion through star-formation may be bypassed. 
However the fact that we do not observe many misaligned systems 
remains. The dynamical timescale (over which gas is expected to re- 
lax into the plane of the galaxy) is shorter for smaller gas reservoirs, 
and thus the observation that Virgo cluster galaxies have smaller 
molecular reservoirs can perhaps somewhat help to explain this dis- 
crepancy. 

Alternatively, if the gas in Virgo galaxies were regenerated 
through stellar mass loss (and thus created kinematically-aligned 
with the stars), then this mass loss would have to come preferen- 
tially from stars at the very centre of these galaxies (perhaps be- 
cause the youngest stars, which lose the most mass, are formed in 
a central starburst as the galaxies enter the cluster). Alternatively, 
the mass loss from the outer parts has to be swept away by ram 
pressure or transported inwards (e.g. by gravitational torques from 
interactions, or by a bar). 

6.2 Molecular gas surface brightness profiles 

In Section [4] we presented the molecular gas radial surface bright- 
ness profiles of the ATLAS''^ ETGs, and compared them with the 
stellar profiles. We showed that, broadly speaking, they can be sep- 
arated into four classes, discs: excess (with or without truncation), 
rings and composite systems. These four classes are very different, 
and hence it is meaningful to ask how these differences arise. 

Molecular reservoirs with regular profiles are almost all found 
in galaxies where the molecular gas is co-rotating (or exactly 
counter-rotating) with respect to the stars (the kinematic misalign- 
ments are tabulated in Paper X). The exceptions are NGC2768, 
NGC7465 and UGC09519, where the molecular gas is in a po- 
lar disc (see |Crocker et al.|2008, A12). HI is detected in 66% of the 
(observed) field galaxies with regular profiles- a higher detection 
rate than observed in the parent population (Paper XIII). This HI 
is often also in relaxed discs (class D or d in Paper XIII), suggest- 
ing many of these galaxies are relaxed systems even out to larger 
radii. This also suggests galaxies with regular profiles have not had 
a large merger or accretion event recently, which would disturb the 
molecular (and atomic) gas disc. 

A third (12/41; 29%) of the galaxies in our sample show 
a CO profile that does not fall off as fast as the stellar pro- 
file. These systems appear to be divided into two types, those 
that show truncations in their molecular gas distribution (7/12; 
58%), and those that fall off regularly but with a different scale- 



length than that of the stars (5/14; 42%). The truncated galax- 
ies either have messy/counter-rotating gas with signs of ongo- 
ing interaction (IC719, NGC 1222; NGC4I50, NGC4550; A12, 
[Beck, Turner & KIoosterman|[2007{ [Crocker et d^[2009) or they 
have evidence of spiral structure or rings in the molecular gas 
(NGC 3182, NGC 3626; A 12, NGC 3489, [Crocker et all[20T0) . 
Similarly, the galaxies which have excess profiles without a clear 
truncation are often messy (IC676, ICI024; AI2, NGC 4753, 
Steiman- Cameron, Kormendy & Durisen|I992) , or may have un- 
resolved rings (NGC4II9; AI2, NGC 4526; Davis et al., 2012, in 
prep). 

The excess galaxies (including those with and without trun- 
cation) are generally found to be low-luminosity ETGs (83% have 
luminosities >-23.5 Afjjg). Lower luminosity (hence presumably 
lower mass) galaxies have shallower potentials, are less able to re- 
tain their gas, and may require longer to force accreted molecular 
gas into discs with similar scalelengths as the stars. All of the ex- 
cess class of galaxies in the field which were observed in HI were 
detected, and many have dynamically and morphologically unset- 
tled gas (as defined in Paper XIII) at large radii. HI at these radii 
takes longer to relax than the centrally concentrated molecular gas, 
and (as argued in Paper XIII and |Oosterloo et ar 2Q10 1 the presence 
of unsettled HI gas suggests recent gas accretion or mergers. Such 
events may thus be the cause of the extended/truncated profiles in 
field galaxies. 

Galaxies with a ring may betray the presence of resonance sys- 
tems, potentially caused by a bar. The molecular gas will then ac- 
cumulate onto stable, non-intersecting closed orbits, as often seen 
in spiral discs (e.g. i Combes| 199 1) . All of these objects show signa- 
tures of bars in their optical morphology or stellar kinematics (Paper 
II). By inspection of the CO PVDs (Section [6. 3^ , we suspect that 
the galaxies that we have classified as having rings based on the 
surface brightness profile generally have gas around the turnover 
radius (often coincident with co-rotation, i.e. inner rings; [ Mazzucaj 
|et al.|20l"T) or further out. Gas within the inner Lindblad resonance 
(i.e. nuclear rings or spirals) is unlikely to be resolved enough to 
show a central dip at our surveys spatial resolution (e.g. the gas in 
NGC 41 19 is likely to be distributed in an nuclear ring given its mor- 
phology, but this is not apparent in its surface brightness profile). It 
has been shown that minor mergers can also potentially cause gas 
to accumulate in rings without the formation of a bar (e.g. |Eliche-| 
[Moral et al.|20TT} . 

The systems with composite profiles are NGC 4710, 
NGC 4694 and NGC 5866. NGC 4694 is a disturbed system with 
the majority of its molecular gas outside the body of the galaxy 
(AI2), and this leads to the derived profile shape. The other 
galaxies in this class have X-shaped PVDs, that usually arise 
from gas caught in two rings within barred galaxy discs (with 
little gas inbetween), the first being gas caught in X2 orbits, and 
the second gas beyond co-rotation (e.g. .Bureau & Athanassoula| 
[1999[ [Athanassoula & Bureau[[1999l i. As discussed in Section |4| 
such profiles are not present in the BIMA-SONG sample, which 
was selected to avoid edge on systems. These galaxies would 
have two clearly separated peaks visible in the surface-brightness 
profiles if viewed at lower inclinations. Thus this types of surface 
brightness profile is not intrinsic to these galaxy, but arises due to 
projection effects. Both systems are edge on, discy fast-rotators, 
and NGC 47 10 has a peanut-shaped bulge (typically thought of 
as being caused by the buckling of a bar; see Bureau & Free man] 
[1999| l, while NGC 5866 is thought to be unbarred with a prototypi- 
cal classical bulge. If the X-shaped PVD of NGC 5866 is indeed 
due to material caught in barred orbits, then this interpretation 
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would have to be revisited. It is possible that we are seeing material 
in inner and outer rings, not associated with a bar, although the 
creation of this double-ring structure may then be difficult to 
explain self-consistently. 

6.3 Velocity profiles 

In Figure [6] (and Appendix B in the online material) we presented 
the CO PVDs of our sample galaxies, overlaid with the modelled 
circular velocity curve, and the stellar and ionised gas rotation 
curves extracted from the SAURON data. We also presented mod- 
eled PVDs, where we use the JAM model circular velocities and 
inclinations, and attempted to fit the observed PVDs with models 
of a relaxed molecular disk or ring (see Section[5^. 

These simple gas surface brightness profiles perform well in 
general, and allow us to test the JAM model circular velocity pre- 
dictions. As revealed in Section |4] however, our real galaxies can 
have quite complicated surface brightness profiles, and thus we do 
not expect to fully reproduce the structures seen in the PVDs. In 
some objects we note that the central parts of the PVD show a 
gap, which is not present in the model data. It seems likely that 
the CO does not extend into the very centre of some of these galax- 
ies (e.g NGC 4435, NGC 4459, NGC 6014). Such central molecular 
deficiencies have been observed in spiral galaxies (e.g. ~ 50% of 
the BIMA-SONG galaxies [Regan et al.|2001^ . In barred galaxies 
these depressions could be due to gravitational torques forcing gas 
outwards to the inner-lindblad resonance, or inward onto a central 
black-hole (for the gas within the dynamical sphere of influence; 
|Combes|200ll . 

6.3.1 Circular velocities 

In Sections [5.1l and l5^ we described the behaviour of the observed 
CO PVDs, and how they compare to the predicted JAM circular 
velocities. The molecular gas reaches beyond the turnover of the 
circular velocity curve in 68±6% of the sample ETGs, as predicted 
from the single-dish spectra in Paper V. This confirms that using the 
shape of the single-dish line profile is a reliable method to select 
galaxies with have molecular gas beyond the turnover radius. 

We separated the galaxies into five rough classes based on the 
relation between the observed CO kinematics and the JAM models: 
those in which the model agreed well with the observations, those 
in which the shape of the model was similar to the observed pro- 
file but predicted higher velocities, those which have polar rings, 
or disturbed gas (and thus we don't expect agreement between the 
model and the observed kinematics), and lastly those in which we 
do not consider due to lack of resolution or known problems with 
the models. 

By comparison with the results of Paper X we find that all the 
galaxies where the models agree well with observed CO profiles 
have their gas kinematics aligned with the stellar kinematics (15/17 
galaxies), or exactly counter-rotating (2/14 galaxies). The CO thus 
rotates in the plane in which the model predictions are valid. Over 
half of these systems have regular CO surface brightness distribu- 
tions (Section|4|. The JAM circular velocity curves trace the outer 
envelope of the CO PVD in many cases. This behaviour is repro- 
duced by our models which include the effects of beam smearing, 
and is often seen in HI PVDs (e.g. Fig 6 in |Oosterloo, Fratemali| 
|& Sancisi|2007 t. In these galaxies with good models the average 
percentage offset between the peak of the observed CO flux at each 
position (the PVD 'trace' introduced in Section [5TT) and the model 
PVD trace is only -3% (~ -5km s~^ ; smaller than the observational 
channel width), with an RMS scatter of 15%. Thus we consider the 
observed CO kinematics and the predicted circular velocity curves 



of these systems to be consistent. The CO in these systems is thus 
likely to be dynamically cold and an excellent tracer of the circular 
velocity. 

In 20% (7/35) of CO-mapped galaxies the shape of the model 
PVD is similar to that of the observed CO rotation curve but is 
at systematically higher velocities. The average offset between the 
CO trace and the model trace is -0.55 dex, larger than would be 
expected from our models which include the effect of beam smear- 
ing. The CO PVDs show that the gas is regularly rotating, and half 
of these systems have regular surface brightness profiles (which we 
argued in Section [6^ are present in relaxed systems). Such large 
differences are thus unlikely to be caused by a dynamical distur- 
bance (e.g. a merger). These systems all have co-rotating or ex- 
actly counter-rotating molecular gas, so the model predictions are 
for the correct plane. These galaxies are also not offset from the CO 
TuUy-Fisher relation for ETGs (presented in Paper V) suggesting 
the gas is relaxed and traces the underlying potential. The remain- 
ing disagreement between the model and observations must be due 
to other factors. Below we discuss possible causes for these remain- 
ing offsets, either due to physical effects within the gas, or potential 
problems with the modeling technique. 

A simple explanation of the offsets between the observations 
and models would be if the molecular gas was not dynamically 
cold. Using the zeroth order approximation shown in equation 9 
reveals you would require very large molecular gas velocity dis- 
persions (aco ~ 80 km s~^) to explain an offset from the JAM 
circular velocity of the magnitude observed in these systems. We 
used the KinMS tool to simulate galaxies with such large velocity 
dispersions, and although the outer-envelopes of the PVD distribu- 
tions can be made to agree, the PVD shape is inconsistent with that 
observed. We also consider such large velocity dispersions to be un- 
likely considering the average val ue for cold atomic gas in nearby 
spirals is only ani ~ 8 km s~^ jvan der Kruit & Shostak||l982 
[Dickey, Hanson & Helou|1990^ . 

Other potential causes of the offset are beam smearing, the 
finite thickness of the gas disc, non-circular motions in the gas (e.g. 
caused by bars) and changes in the density profile (e.g. spiral arms), 
which could all act to lower the gas (apparent) rotation velocity. The 
first two effects are included in our models, and are thus unlikely to 
be the cause of the offsets, and thus we concentrate on the other 
possible causes here. 

Around half of the systems in which the JAM model circular 
velocities are higher than the gas velocity have central CO rings 
(NGC 3182), or have (single or multiple) solid-body components 
in their observed PVDs (NGC 2764, NGC 47 10, NGC 5379). Such 
morphological and kinematic features often betray the presence of 
resonance systems, potentially caused by a bar (e.g. [Combes] 1991) . 
The majority of the galaxies (6/7) also have a bar visible in opti- 
cal imaging (Paper II). In barred systems material follows elliptical 
rather than circular orbits, and hence some deviation from the cir- 
cular velocity is expected. Non-circular motions of this type can 
provide a possible explanation for the offset between the observa- 
tions and models in galaxies with bars. One must also be cautious, 
as material caught in a circular ring at a resonance point such as 
co-rotation does rotate at the circular velocity. Without understand- 
ing the resonance structure of each galaxy it is thus difficult to say 
with certainty which galaxies have lower CO velocities due to non- 
circular motions in the gas. In principle one can detect non-circular 
motions by the effect they have on the observed molecular gas ve- 
locity fields (see e.g. [Schoenmakers, Franx & de Zeeuw|1997[[Kra-| 
[jnovic et al.[2006^ , however our modest spatial resolution limits our 
ability to do this with the current data set. 
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An example of a galaxy which may have gas in non-circular 
orbits is IC 0676. This galaxy has a strong optical bar, and the CO 
rotation is aligned with the bar, not the main disk. As the PVD of 
IC 676 does not show solid body rotation it is unlikely to have a 
molecular ring, and hence the gas velocity could be affected by non- 
circular motions, and produce the observed discrepancy. 

Lablanche et al., 2012 (hearafter Paper XII) have also shown 
that the presence of a bar can induce an additional uncertainty (of 
the order of ^15%) in the M/L predicted by the lAM model. 
As the assumed M/L sets the normalisation of the rotation curve 
{v oc M/L), this could explain some (but likely not all; the av- 
erage effect would be ~15 - 20 km s~^) of the offsets of the lAM 
predictions. As we have 7 galaxies with JAM profiles that predict 
higher velocities than the CO, and few where the predicted circular 
velocity is low, statistically it is unlikely that this is the dominant 
cause of the difference between the observed CO rotation and JAM 
predictions. 

The models we show here are based on the JAM circular ve- 
locity curve for an edge on system. These were derived from mass 
models of the galaxies, and the observed SAURON kinematics. The 
output circular velocity curves have been projected back into the 
line of sight with the best fitting JAM inclination, in order to allow 
comparison with the observed CO PVD. If the JAM models system- 
atically underestimated the inclination of the system then this could 
result in the model predictions being too high. Careful comparison 
with dust distributions in these systems, and other inclination indi- 
cators (e.g. those presented in Paper V) suggests that the derived 
JAM inclinations are unlikely to be biased in this manner however 
(Cappellari et al., in prep). 

Although the early-type galaxies in this sample were selected 
to show no evidence of spiral arms in optical imaging, it is pos- 
sible for them to have spiral structures in the molecular gas disk 
(c.f. NGC 3489; [Crocker et al.|2011^ . Such structures may remain 
unnoticed if the beam size of our observations is too large to re- 
solve them, or if the gas disk is observed close to edge on (5/7 of 
the galaxies where the observations and models do not agree have 
inclinations >60°; Paper V). Both spiral structure itself, and asso- 
ciated streaming motions along spiral arms could have a noticeable 
effect on the observed kinematics of the gas, and the PVD. Stream- 
ing motions can have a large amplitude (e.g. ~50 km s^^in M51; 
[Shetty et al.|2007| , however these systems would have weaker spi- 
ral structure than M51, and at the scales probed here the measured 
velocity should change by <10km s~^, not enough to explain the 
offsets we see. 

The effect of spiral structure itself (e.g. not having gas at all 
radii along the PVD) depends on how tightly wound the spiral arms 
are, and the angle they make to the line of sight. Figure [TO] shows a 
simulation (made in the same fashion as described in Section [5^ 
of the PVD observed for three different spiral gas distributions. The 
first column shows a simple disk simulation as a reference, while 
the second column shows a tightly wound spiral, and the third and 
fourth panels show a loose spiral, viewed at different angles. For 
each simulation we show the face on view of the input surface 
brightness profile in the top row, and how it would look when pro- 
jected to an edge on inclination in the second row. The bottom row 
shows the observed PVD for the edge on system in black, overlaid 
on the PVD of the disk system. A tightly wound spiral pattern dis- 
plays a very similar PVD to a complete disk, and hence it seems 
unlikely that a tightly wound spiral could cause the offset between 
simulations and observations that we observe. A loose spiral pat- 
tern however displays a PVD that changes in form depending upon 
the angle at which it is viewed. It is important to note however that 



the observed PVD in these cases covers some fraction of the to- 
tal area of the PVD of the disk simulation, but does not go outside 
it. One can thus potentially identify systems with such spirals by 
comparison with the PVD of a disk (as presented in Figure |6] and 
Appendix B in the online material). None of the cases in which 
the JAM and CO rotation disagree are clear candidates for having 
such spiral structure, based on the PVD, but we are unable to rule 
out such an explanation completely. Even if some objects do have 
such a structure, the amplitude of this effect would not be enough 
to explain all the offsets seen. 

A potential cause for the models circular velocities to not agree 
with the observed CO rotation is the effect of dust on the MGE. 
Many of the systems where the JAM models overpredict the CO 
velocity have central dust obscuration in the optical images used to 
create the MGEs (even after a dust correction based on g-r colour 
is applied). Such system include IC 676, NGC 2764 and NGC 4710. 
The models of these galaxies are likely to badly capture the light 
distribution in the very central parts, that are vital to set the circular 
velocity normalisation. The measured light profiles may become 
too shallow in the inner parts and the models will require a high 
M /L to best fit the (steeper) rise in the stellar velocities. 

Assuming a fixed dust-to-molecular gas ratio, it is clear that 
we would expect these problems to be worse at high molecu- 
lar gas masses. Systems with profiles that are at higher veloci- 
ties than the observed CO rotation do tend to have larger molec- 
ular gas reservoirs (Paper IV). These galaxies have a mean H2 
mass of 6x10*Mq, nearly half an order of magnitude higher 
than that of galaxies with good JAM model predictions (< 
Mh2 >~2x IO^Mq). Such problems with dust could be overcome 
to some degree by utilising near-infrared imaging (which is less af- 
fected by dust than optical images) in the centres of these objects 
when creating the MGE, but creating such models is beyond the 
scope of this paper. Comparison of the rotation curve of NGC 47 10 
derived by [Williams, Bureau & Cappellari| |2009| using infrared 
imaging shows that an offset of «30 km/s can be induced by this 
effect, enough to explain the offset between the JAM velocity and 
CO in this object. 

One other potential cause of the offset between the models and 
observations would be if the assumption of a constant AI/L as a 
function of radius is invalid. The systems with model PVDs that lie 
above the CO profiles tend to have large H/3 emission equivalent 
widths, and low [0III]/H/3 ratios (see Figs.|7]and|9j. These systems 
also have larger molecular masses, on average, as stated above. The 
EW of H/3 is often used as a star-formation tracer in spiral galaxies 
(where an EW(H/3) <; 2A is usually considered indicative of on- 
going star-formation). One must be cautious when using EW(H/3) 
to trace star-formation in ETGs, as other sources of ionisation can 
cause H/3 emission (e.g. old stellar populations, [Sarzi et al.|2010) . 
However, in all but one of the galaxies here that have EW(H/3) 
_^ 2A the ionisation is likely to be from star-fo rmation (as Figure 



Sarzi et al. 



2010 1. 



I^shows they have log([OIII]/H^) ratios < -0.2; 
The one galaxy where this is not the case is NGC 7465, which is a 
composite source, with ionisation both from star-formation, and a 
X-ray bright AGN. 

Active star-formation would create a population gradient, with 
many young stars in the inner parts of the galaxy where the CO 
is concentrated, and older stars dominating in the gas-poor regions 
further out. The galaxy average M/L (and thus circular velocity) 
predicted by the JAM model would thus be too high in the inner 
parts of the galaxy where the CO is, leading to an overestimation 
of the circular velocity (with a corresponding underestimation in 
the outer parts), [^ung. Bureau & Cappellari[ ( [2008^ show that the 
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magnitude of this effect can be as high as a factor of 3 in velocity in 
the inner parts of NGC 3032, a galaxy with a strong population gra- 
dient. Even in galaxies where the population gradient is shallower 
this effect could be large enough to explain the difference between 
the observed and predicted circular velocities. Cappellari et al., in 
prep identify galaxies with H/3 absorption linestrengths > 2.3 A 
(within an effective radius) as likely to have population gradients 
strong enough to affect the JAM models. Of the systems classified 
as having JAM model circular velocity profiles that predict higher 
velocities than the observed CO rotation, 5/7 have such strong in- 
dicators of young stars (IC0676, IC0719, NGC 2764, NGC 5379, 
NGC 6014). The two galaxies with JAM models which match well 
the molecular gas PVD but which have large H/3 emission equiva- 
lent widths are IC1024 and UGC 06719. In these galaxies the CO 
extends over the entire disc of the galaxy, to beyond 1 Rc, and hence 
the variation of M/ L as a function of radius is likely to be lower. 

Figure [8] shows that many of the systems in which the JAM 
model and CO rotation do not agree are at low velocity dispersions. 
If the difference between the JAM model velocities and the CO rota- 
tion is due to ongoing active star-formation then this could be a con- 
sequence of 'downsizing', where lower mass galaxies have more 
star-formation activity at low redshifts (e.g. |Cowie et al.]|1996^ . 
Even if the star-formation rate of all ATLAS'''^ galaxies were the 
same, the luminosity weighted M/ L in a low-mass galaxy would be 
more affected than in a high-mass galaxy. As we touched on above, 
however, the molecular gas reservoirs of these low-mass systems 
are much larger on average, so the situation is even more extreme. 

We thus conclude that in the galaxies with the largest molecu- 
lar masses, and highest star-formation rates, dust obscuration, non 
circular motions in the gas and/or population gradients could be af- 
fecting the observed CO, or derived JAM circular velocities. The 
magnitude of these effects can be up to 0.6 dex in velocity. In the 
context of the entire ATLAS"^'^ sample however, < 3% of galaxies 
(8/260) have sufficient molecular gas, dust and star formation to af- 
fect the JAM model predictions, so these effects should not strongly 
bias results for the whole sample. Indeed, all cases without the prob- 
lems highlighted above we found excellent agreement between the 
CO observations and JAM models predictions. 

In the 11% (4/35) of CO-mapped galaxies where the CO is 
disturbed (as identified in AI2), its rotation (or lack thereof) is un- 
related to the JAM prediction of the circular velocity. In one of 
these systems (NGC 1 222) the JAM model is also likely to be bad 
due to large scale stellar kinematic disturbances, however in the 
other cases the stars seem unaffected by the gas disturbance. All 
of these galaxies have extended surface brightness profiles, with 
truncations, supporting our conclusion that these profiles indicate 
mergers and/or recent accretion. 

Three of the CO-mapped galaxies (9%) have their molecular 
gas rotating in the polar plane (Paper X), hence the disagreement 
with the model predictions (for rotation along the photometric ma- 
jor axis of the galaxy) probably implies that the total mass distribu- 
tion in these three systems is flattened (as a spherically symmet- 
ric mass distribution would have the same rotation curve in any 
plane). The luminous matter in these three galaxies has elliptici- 
ties between 0.45 and 0.65, supporting this assumption. If one has a 
suitable tracer (such as HI) that extends to large enough radii it may 
be possible to constrain the shape of the dark matter halo in these 
systems (e.g. as done for IC2006 in |Franx, van Gorkom & de Zeeuwj 
|1994^ , however with this data such an analysis is impossible. 

6.3.2 The dark matter halo of NGC4753 

|Steiman-Cameron, Kormendy & Durisen| ( |I992| > have discussed in 



detail the complex, twisted dust lanes of NGC 4753, and showed 
that they are consistent with a disc that is strongly twisted by differ- 
ential precession. They used this twisted disc to set constraints on 
the shape of the matter distribution, finding that the total dark and 
luminous mass is nearly spherically-distributed, with an ellipticity 
no higher than 0.16 and no smaller than O.OI. Due to a lack of kine- 
matic information for the warped disc however, they were unable 
to determine if the dark matter halo is oblate or prolate. If the ma- 
terial in the warped disc rotates in the same direction as the stars, 
then the dark matter must have an oblate shape, while if it counter- 
rotates the halo is likely to be prolate. Assuming that the molecular 
gas we detect in this galaxy is associated with the warped dust disc, 
it is possible to break this degeneracy. The molecular gas is kine- 
matically aligned with the stars and thus it follows, according to 
|Steiman-Cameron, Kormendy & Dur isen ( 1992 1, that the dark halo 
in this galaxy is oblate, with a flattening between 0.16 and 0.01. 

6.3.3 Ionised gas velocities 

The ionised gas in the CO-mapped ATLAS'^° galaxies always 
(within the errors) rotates at the same speed, or slower, than the 
molecular gas. This confirms that CO is the best tracer of the circu- 
lar velocity in the inner parts of ETGs. Due to the large mismatch 
in angular resolutions between the CO interferometry and the HI 
data we possess (Morganti et aL|2006||Oosterloo et al._20I0 , Paper 
XIII) we are unable to directly compare CO and HI velocity pro- 
files in this work, but we would expect HI to be a comparable tracer 
of the circular velocity, and the best available probe at large radii. 
Comparisons of the integrated profile shapes will be presented in a 
future work, discussing the HI TuUy-Fisher relation of ETGs. 

The differences between the velocities of the CO and ionised 
gas relates to the ionisation mechanism powering the optical emis- 
sion. HII regions embedded in a dynamically cold star-forming 
molecular disc are likely to rotate much like the disc itself, and be 
dynamically cold. On the other hand, gas clouds ionised by other 
(most likely diffuse and stellar; [Sarzi et al.|2010^ , sources do not 
have to necessarily display such a relaxed kinematic behaviour. 
In the specific case of emission powered by shocks between gas 
clouds, for instance, a larger velocity dispersion for the ionised- 
gas is expected. In Figure |9] it is clear that the most star-forming 
galaxies (with the largest H/3 equivalent widths and star-forming 
[0III]/H/3 ratios; Figure |9]( have the smallest differences between 
the ionised gas and the molecular gas velocities. 

A simple first order estimate based on equating the second mo- 
ments of molecular and ionised gas (as in Equations 8 & 9) shows 
that the small systematic offset from the dashed zero- velocity offset 
line of Figure |9] can be explained by the observed velocity disper- 
sion providing dynamical pressure-support to the ionised gas (A 
AV ~21 km s^^is expected given a < agas > ~ 62 km s^^). 
Galaxies with a smaller amount of star formation have a wide range 
of offset values between the ionised and molecular gas velocities, 
reflecting a variety of ionisation mechanisms and the presence of 
few, if any, gas clouds in the immediate vicinities of OB-stars and 
their giant molecular nurseries. 

7 Conclusions 

We have investigated the extent and surface brightness pro- 
files of the molecular gas in the CO-rich ATLAS^° ETGs, and 
compared them to spiral galaxies. We find that the molecular gas 
extent is smaller in absolute terms in early-type than in late-type 
galaxies, but that the size distributions appear to scale similary with 
optical/stellar characteristic scalelengths. Amongst ETGs, we find 
that the extent of the molecular gas is independent of the kinematic 
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Figure 10. Simulated observations showing the effect of spiral structure on the derived position-velocity diagrams. The left column shows a disk simulation, 
and its PVD for reference. The second column shows a tightly wrapped spiral like the Milky Way (with a pitch angle of 12°; |Vallee|1995) , and the third and 
fourth show a loosely wrapped spiral patterns (pitch angle 35°) observed with different angles to the line of sight. The top row shows the input surface density 
distribution in each simulation, viewed face on. The same radial suiface brightness profile is used in each case. The middle panel shows the zeroth moment of 
the simulated observation, viewed edge on. The bottom panel shows the disk PVD (colours), and in the last three panels this is overlaid in black with contours 
of the PVD for the spiral density distribution. The blue line shows the input circular velocity curve, which is identical in each case. The simulations all were 
performed with a resolution of 4'.'5, a pixel size of l"x l"and a gas velocity dispersion of 8 km s~^. 



misalignment, despite the many reasons why misaligned gas might 
have a smaller extent. The extent of the molecular gas does de- 
pend on environment, with Virgo cluster ETGs having less extended 
molecular gas reservoirs. This may be due to ram pressure stripping, 
gravitational interactions within the cluster and/or a greater preva- 
lence of bars driving gas towards the galaxy centres. 

We have identified four classes of surface brightness profiles 
in the CO-mapped ATLAS'^° ETGs. Around half have molecular 
gas that follows the stellar light profile, similarly to molecular gas 
in spirals. These galaxies appear to have relaxed gas out to large 
radii, suggesting that these galaxies have been left undisturbed long 
enough for the molecular gas to relax into a profile that follows 
the stellar mass. A third of the sample galaxies show molecular 
gas surface brightness profiles that fall off slower than the light, 
and sometimes show a truncation. These galaxies often have a low 
mass, and either have disturbed molecular gas, suggesting that a re- 
cent merger or gas accretion may cause the excess, or are in cluster 
environments where ram pressure stripping and the presence of hot 
gas may compress and/or truncate the gas. The remaining galaxies 
have rings, or composite profiles, that we argue can be caused by 
the effects of bars (or potentially minor mergers). 

We also investigated the kinematics of the molecular gas using 
position-velocity diagrams. On to these we overlaid the JAM model 
circular velocity predictions and the stellar and ionised gas veloc- 
ities. We find that the molecular gas reaches beyond the turnover 
of the circular velocity curve in 75±6% of our ATLAS^^ ETGs, 
validating previous work. The CO is also closer to the predicted 



circular velocity than the stars or ionised gas, validating its use as 
a kinematic tracer for TuUy-Fisher and similar analyses. This study 
suggests that the dust lanes in NGC 4753 are likely to co-rotate with 
the stars. From this it follows, according to |Steiman-Cameron, Ko-| 
[rmendy & Duri sen (19921, that the shape of the dark-matter halo in 
NGC 4753 is oblate, with a flattening between 0.16 and 0.01. 

We found that in ~50% of our galaxies the predicted circular 
velocity profiles and the observed CO rotation are consistent, the 
CO appears to be dynamically cold, and CO is an excellent tracer 
of the circular velocity of the system. In 20% of CO-rich galaxies 
we found the JAM model prediction of the circular velocity is sys- 
tematically above the observed CO velocity. These galaxies have 
strong dust features, that make constructing a mass model from op- 
tical imaging challenging. These systems also often have bars, and 
higher molecular gas masses and H/3 equivalent widths than the 
rest of the sample. This suggests that dust obscuration and/or pop- 
ulation gradients can affect model predictions of the circular veloc- 
ity. At the same time bars, or non axisymmetric structures could 
cause the gas to rotate slower than the circular velocity (however 
the amplitude of this effect is likely insufficient to explain the ob- 
served discrepancies). As systems with CO that show the problems 
highlighted above make up a small percentage («3%; 8/260) of the 
whole ATLAS''^ sample this will not strongly bias results based on 
JAM models of this type when applied to all galaxies. In all cases 
where the CO appears to be regularly rotating, where the MGE 
models well reproduce the observed photometry, and where star for- 
mation activity is low we found excellent agreement between these 
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two completely independent measures of the circular velocity (CO 
observations and JAM models of the SAURON stellar kinematics), 
in agreement with an earlier study on a smaller number of objects 
by |Young, Bureau & Cappellar i (2008). 

We also note that the ionised gas in the most star-forming 
galaxies appears to be dynamically colder, likely because it is as- 
sociated with star formation in the molecular disc. 

In conclusion, we find that the molecular gas reservoirs in 
many ETGs are very similar in terms of their relaxed kinematics 
and relative extent to those in spiral galaxies. However rich environ- 
ments and merger activity can lead to some systems having different 
morphological and kinematic properties. Understanding how these 
differences affect the physical conditions within the gas, and thus 
star-formation, will be required to correctly understand the ongoing 
evolution of the most massive galaxies in our universe. 
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Figure Al. Radial surface brightness profiles of tlie stars (red dashed line; r-band) and molecular gas (blue solid line with errors) for our 
ATLAS'*'^ ETGs (normalized at the first CO datapoint). The CO beamsize is shown as a black dotted line. The error bars on the CO measure- 
ments denote the width of the elliptical annuli in the x-direction and the RMS noise in the elliptical annulus in the y-direction. The letters at the 
bottom-right corner of each plot denote the profile class. "Reg" denoting a regular profile, "E" an extended profile, "E(T)" an extended profile 
with truncation, "R" a ring, "C" a composite, and "U" an unresolved profile. 
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Figure Al. continued 
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Figure Bl. CO Position- Velocity diagrams (orange witli black contours) for the ATLAS^'^ galaxies, overlaid with the JAM circular velocity 
curve (thick black line). On the left, these have also been overlaid with the observed stellar (red dashed Hne) and ionised gas (blue line) rotation 
curve. On the right we overlay the modeled observations of each galaxy projected into the line of sight (blue contours), created assuming that 
the gas has an intrinsic velocity dispersion of 8 km s^^, is rotating at the predicted circular velocity, and has an exponential or gaussian surface 
brightness profile, as described in Section 5.2. The MGE model of IC0676 is effected by dust (Scott et al., in prep), and as discussed in the 
text the gas may have significant non circular motions due to a bar. These two effects likely cause the discrepancy. NGC 0524 has a degenerate 
inclination solution which causes the disagreement between model and observations (see figure Al in Cappellari et al., 2006) 
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Figure Bl. continued. NGC2764 lias an MGE model which is heavily affected by dust (Scott et al., in prep). 
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Figure Bl. continued 
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Figure Bl. continued. NGC4710 has an MGE model heavily affected by dust, and bar, which together may cause the disagreement between 
model and observations (Scott et al., in prep). 
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